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ABSTRACT*  This  is  a survey  of  recent  developments  about 
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istics  are  discussed  for  such  a flow  in  prismatic  ducts  with 
symmetric  cross  sections  (closed  or  open;  as  well  as  in  the 
boundary  layer  of  a flat  plate. 


Hovr  does  the  ratio  ( (*  • ) of  turbulent  to  molecular 

shear  stress  vary  with  a dimensionless  distance  ( vi  » /ay*/*) 
from  the  wall?  Of  late,  in  reply  to  this  question,  the 
function  /*fa)  has  been  approximated  analytically  (References 
(m) , (q)  ,f  (tn  over  several  regions  of  the  boundary  layer. 
This  allowed  to  trace  as  well,  by  integration,  a wean  value 
profile  respectively  for  the  velocity,  and  the  temper a ture. 
directly  as  to  the  former,  step-wise  for  the  latter. 


In  the  first  part  of  this  report,  physical  requirements  of  the 
lower  atmosphere  are  transposed  into  analytic  boundary  conditions 
for  the  mean  velocity  profile.  It  is  discussed  in  regard  to 
a current  formulation  with  two  constants  - power  relation, 
logarithmic  law  - anci  also  in  respect  to  more  recent  endeavours 
(H.  Reichardt,  K.  Elsser)  for  an  analytic  expression  of  the 
continuity  over  the  whole  border  film. 

The  early  dichotomy  (Prandtl,  Taylor,  etc.)  of  the  boundary 
layer  evolved  through  a subsequent  trichotomy  (von  Kdrmdn, 
iteicnardt,  etc*)  to  the  latest  theory  (Reichardt,  Eiser,  etc.). 
This  is  reflecting  e continual  decay  of  turbulence  effects 
from  the  free  stream  to  their  extinction  at  the  wall. 

Reichardt *s  unified  velocity  profile  (1951)  is  set  forth  and 
discussed  ns  to  the  prevalent  term,  near  the  wall  and  in  a 
median  region  respectively,  and  as  to  combining  the  terms  in  a 
complete  profile." 
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profile  admits  of  a rigorous  computation* 
field  of  turbulence  no  general  analytical  solution  is  known. 
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Yet,  for  the  thermal 


The  second  part  deals  with  progress  toward  solving  the  turbulent 
case  by  integrating  the  temperature  equation  numerically  or 
graphically  .in  successive  steps.  One  or  two  steps  appear 
sufficient  for  air  and  fluids  with  higher  Prandti-nurabers 
(cils,  etc.),  while,  at  most,  four  steps  would  be  needed  in 
ease  of  very  low  Pr-numbera  (molten  metals  or  alloys  etc.#)* 

The  third  part  deals  mainly  with  an  iterative  integration  as 
a meana  of  solving  the  heat  transfer  problem  for  turbulent 
flow.  Two  integrating  varieties  of  ths  analytical  method  are 
discussed. 

Theoretical  results  are  compared  with  those  of  global,  semi« 
empirical  relations  for  the  heat  transfer  coefficient. 

The  investigated  results  cover  p.  wide  range  of  the  generalized 
Prandtl-number ; 0 ^ ^ jo* 

Numerical  values  of  the  heat  transfer  coefficient,  as  obtained 
by  successive  integration,  appear,  currently,  in  good  agreement 
with  those  based  on  semi-empi rical , product  relations  of 
po7/ered  similarity  parameters  as  factors  ( Reference o : ( f ^ ), 

von  Karman,  and  (ra*  ),  Stei  chardt*  for  high  Pr-numbers; 

(q  ),  Riser , ana  Is),  Lyon,  for  low  Pr -numbers ) . 

The  fourth  part  Is  presenting  soma  concluding  remarks. 
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Tbi  a NAVOPD  is  a review  of  recent  progress  on  the  influence 
or  turbulence  upon  heat  transier  at  3mooth  walls  for  a flow 
in  prismatic  ducta  or  along  a flat  plate.  Suet  developments 
concern  problems  currently  encountered  in  the  NCL  AeroBailiatic 
Research  D«partment.  The  more  recent  of  this  reviewed  rec  jarch 
work  was  published  in  Gottingen  and  ZUrioh.  No  translatioh  in 
English  is  eveilablo  at  present.  Therefore,  it  is  felt  worth 
while  to  include  and  correlate  the  new  methods  and  results  in  a 
comprehensive  presentation. 

This  investigation  was  carried  out  under  task  number  Re9a-108-1. 


EDWARD  L.  WOODTARD 
Captain,  USN 
Comma  nder 


H.  H.  KURZWEG,  Chief 
AeroballiiStic  Research  Department 
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Symbols 

Coordinates  (x  parallel  to  wall*  y normal 
to  wall,,  z normal  to  x-y  plane) 


Mean-temporal  components  of  velocity  In 
direction  of  coordinates  x,  y,  z 


‘rVl 


Inst an t an ecus  values  of 
of  velocity  components 


d ® rr 

A *•*!/  - — " 


Instantaneous  values  of  resultant  velocity 
components 


Root  mean  square  vr  m s)  values  (or 
effective  values)  of  the  fluctuating  parts 
of  velocity  components.  For  fluctuations 
of  sufficiently  high  frequency,  & finite 
time  interval  "tn  can  be  large  enough  for 
establishing  quadruti ; mean  values  of  the 
fluctuation,  while  the  mean  velocity  may 
be  assumed  as  practically  constant. 


A • 


Cross-sectional  area. 


ajl  » A \ziiA 

v s 


Average  value  of  mean  velocity  over  cross- 

^4  -.1  

w Awuax  wroti* 


U 


Largest  value  of  u-oomponent  of  velocity. 


T Absolute  temperature. 
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T-Ta  Difference  between  temperature  of  moving 

fluid  and  wall  temperature. 


0 


Largest  temperature  difference. 


£[■  Dimensionless  temperature  difference. 

T - A*'  l — Mean  temperature  of  moving  fluid  over 

**•  J/-  cross  section. 


f 


Density  of  roving  fluid 


r 


js  u .9  • Dynamic^ and  kinematic  coefficient  of 

molecular  viscosity. 


Coefficient  of  apparent  viscosity  due  to 
turbulence . 


A Coefficient  of  heat  transport  by  ouruootiva 

<j,  through  turbulence. 


Ratio  between  coefficients  of  turbulent 
heat  flux  density  and  apparent  viscosity. 


• Coefficient  of  molecular  heat  conductivity. 


j c^aT  Enthalpy. 
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Specific  heat  at  constant  pressure. 


Shear  stress. 


i 


Heat  flow  density  (heat  per  unit  time  and 
r unit  area' . 


Subscript  Rt"  refers  to  turbulent  contribution,  subscript  "m'’ to 
molecular  contribution. 


Pr  and  1 1 = numb  s r 


Generalized  Prandtl-number . It  is  a 
characteristic  parameter  of  turbulent 

< 4*  Vi  V.  j_i 

* n a w**  uw  o w VS  ou iOJ.  U1  « 


T 


st  A4 


U" 


Dimensionless  velocity  referred  to  free- 

a ^rGciu  Vfill)  A * 


Shear  stress  velocity  at  wall. 
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vi,  W * 
^ * l -u. 


Dimensionless  local  velocity  referred  to 
shear  velocity  at  wall. 


Dimensionless  distance  from  wall. 


Dimensionless  coordinate  in  z -directi on. 


*7  «(Vl. 

vyr 


Dimensionless  width  of  viscous  boundary 
defined  as  that  distance  from  the  wall 
beyond  which  turbulence  accounts  for  a 
higher  apparent  shear  stress  than  the  ps 
due  to  molecular  agitation. 


• • • . Dimensionless  width  of  the  heat  conducting 
4 yi.Vc  boundary  lay*r  defined  as  that  distance 
* from  the  wall  beyond  which  turbulence 

accounts  for  a higher  transport  of  heat 
than  the  molecular  conductivity  doea. 

*■  J r,  i xftitM  R®fch*rdt  integrals  concerning  the  temper- 
n r*  v r*-J  a^ura  profile  of  the  turbulent,  heat 


transmitting  boundary  layer, 

S * - .t/jg-  - \ 

1 -rU 


t.  * 


fV 

\ Jk  «•  ' 

I l 


* * ' 


^ » o 


A.  v. 


Radius  of  circular  tube  or  distance  between 
wall  and  symmetry  plane  for  a prismatic 
channel. 
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Dimensionless  either  radius  or  distance  of 
wall  from  symmetry  plane. 


H - i 


Distance  from  axis  of  circular  tube  or  from 
symmetry  plane  for  a prismatic  channel. 


Reynolds  number  referred  to  the  mean 
velocity  over  cross-section  and  to  the  tube 
diameter  or  largest  channel  dimension  across 
the  symmetry  plane. 




UT„-T.) 


Nusselt  number  referred  to  heat  flow  through 
the  wall  surface  and  to  the  average  mean 
temperature  of  tho  moving  fluid  over  a 
cross-section. 


_ 

J r ~r  T'  \ 
K V.T*-  *.) 


• Nusselt  number  referred  to  the  difference 
between  the  equilibrium  temperature  at  an 
insulated  wall  and  the  temperature  of  the 
actual  wall . 


’Th-T. 


Diuieiisiouel  heat  transfer  coefficient  based 
on  the  temperature  difference:  T — T 


* — . . Equilibrium  temperature  at  a perfectly 

c insulated  wall. 


v,  = 


■ Dimensionless  heat  transfer  coefficient 
bt.sed  on  largest  velocity  ( U ) and  largest 
difference  of  temperature  ( Q ). 
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, . . . . Additional  velocity  for  the  median  region 

of  flow  as  introduced  in  new  theory 
(Reichardt,  I95O-.51)  in  order  to  account 
for  a slight  variance  between  experimental 
results  and  the  basic  logarithmic  law. 


Kf  4i , Constant  magnitudes. 

k - 4r  » *. 
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A RFVIEW  OF  INVESTIGATIONS  ON  TURBULENT  FLOW 
WITH  HEAT  TRANSFER  AT  SMOOTH  WALLS 


Prefatory  Remarks 

1.  The  characteristics  of  turbulent  flow  in  prismatic  ducts 
with  cross-sectional  symmetry  or  along  a flat  plate  are  briefly 
surveyed.  This  is  done  with  the  aim  of  a closer  approach  to 
the  heat  transfer  conditions  at  smooth  walls.  The  wall  tempera- 
ture is  assumed  constant  over  the  entire  area  of  the  solid 
boundary. 

2.  The  duct  may  be  either  of  closed  cross  section  (tuba)  or 
open  ( channel) . 

3.  The  fields  of  mean  velocity  and  mean  enthalpy  (or  tempera- 
ture) are  based  on  a more  complota  account  of  the  continuity 
of  flow  and  the  other  physical  requirements  at  the  wail 
boundary  and  the  free-stream  portion  of  ee.ch  profile. 

[}..  The  flow  conditions  are  limited  as  a rule  to  negligible 
effects  on  the  density  of  the  fluid. 
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Part  I«  The  Profile  of  Mean  gelooity 
Introduction 


5.  In  thla  first  part  of  tha  Report,  only  the  velocity  field 
is  considered.  After  a brief  historic  survey,  a discussion 
will  follow  on  Heiehardt**  recent  (1950-51;  achievement  in 
unifying  the  analytic  expression  for  the  complete  velocity- 
prof  lie. 


Boundary  Conditions 


6.  A parallel  mean  flc??  in  redirection  is  considered.  For 
the  profile  of  mean  velocity  in  the  lower  atmosphere,  the 
physical  requirements  at  a smooth  wall  are t 


(*) 

Hon-slip  condition,  i.e.  zero  value  for  the  mean 

velocity:  

AA.  - 0 
o 

( ) 

Finite  shear  stress,  i.e.  finite  gradient;  of  mean 
velocity  i 

( da  \ ° 

\&7  L ^ 90 

it) 

Continuity  equation  extended  to  include  the  three* 
group  of  the  fluctuating  parts  (u.  , ?,  , 
w ) of  the  velocity  components*  for  them:  i £ 

» v 

r1  mi 


r\  - 

V f 


Vs, 
- i 


C-'U// 

Ti 


7*  At  the  symmetry  plane  of  a prismatlcal  chauiel  or  at  the 
axis  of  a circular  tube,  the  physical  requirements  for  the 
velocity  profile  *»v  ^ ; 


<*> 

U> 


I Js 


Vanishing  velocity  gradient: 


1 'c^tVw 


Finite  curvature  ( D at  the  transit ion  point  for 
two  symmetrical  branches  of  the  velocity  profile, 
i.e.  an  analytical  maximum  for  the  velocity  at 
center : 


1 4V 

\ ir° 

V dy* 

Ctk(U 

* A. 

• 

N 


I 


mdi* 
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Power  Relation  . Logarithmic  Law 


6.  It  is  known  that  for  the  flow  near  a wall,  though  not  too 
closa  to  it,  the  measured  velocity  profile  may  be  fipproxim^tod 
by  i 


( a 


A power  relation  with  a fractional  exponent: 
dependent  on  the  Reynolds  number: 


c {•*.!  *h 

v ) i 


<1) 


In  a range  of  the  Blasius  values: 

n = 7,  c •-=  8.7  are  valid  up  to  about  Re  = 100,000. 
Above  this  limit,  the  actual  experiments  account 
for  Increasing  n-values  with  Re,  via.  n — » 8,  9,  10. 


( b ) A logarithmic  law  with  two  constants: 


t>o  + (2) 

L M Z. S’  ( Tcw-w-HC 

9 * \ 
u.  Kfc*-***^  ) 

V 

Tills  formula  would  asea  & more  general  approximation  than 
the  power  law. 

9»  It  la  easy  to  check  that  the  power  relation  ( a)  - though 
giving  non-slip  at  the  wall  - violates  there  the  boundary 
condition  of  a finite  velocity  gradient. 


*1  A Am  a.Viwm  AU*  1 » i.v. m _ a • — j»  . . ...  / *_ 

JLU*  on  tns  ctncr  nana,  tns  lu^u’nimia  x*w,  m 1 ora  v o 
violates  at  the  wall  the  non-slip  condition  and  that  of 
finite  velocity  gradient. 


> 


11.  Besides,  both  laws  (a)  and  (b)  do  not  satisfy  the  center 
conditions  for  a symmetric  duct  or  the  circular  tubs. 


12.  In  order  to  remedy  tha  Inconsistences  in  respect  tc  the 
boundary  conditions  at  the  wall.  It  had  become  customary  to 
assume  that  the  turbulent  agitation  would  not  reech  the  fluid 
that  moves  in  immediate  vicinity  of  the  wall.  Thus,  a purely 
laminar  sublayer  of  very  small  thickness  was  postulated  in  a 
position  between  the  wall  and  the  turbulent  region.  Within  the 
sublayer  a linear  velocity  distribution  vras  currently  assumed, 

ic«;  : 


1 M \ 


■n 


(3) 
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13.  This  implied  a constant  velocity  gradient,  and  hence  an 
invariable  shear  stress  from  the  wall  through  the  thickness  of 
the  sublayer.  Yet,  as  a matcsr  of  experimental  evidence,  the 
shear  stress: 


diminishes  with  increasing  distance  from  the  wall.  To  account 
for  this  fact,  formula  (3)  ooulu  of  course  be  replaced  by  a 
parabolic  branch  (Poir.eu.ille)  to  hold  in  the  sublayer. 

lli , There  are  cogent  considerations,  however,  that  seem  to 
preclude  the  very  existence  of  a purely  laminar  sublayer* 

For  a flow  along  a solid  wall,  it  seems  physically  inconceivable 
to  limit  the  spreading  of  turbulenea  to  an  arbitrary  fluid 
boundary  of  a purely  laminar  sublayer.  It  would  appear 
questionable  as  well  to  associate  the  fully  developed  turbulence 
to  a local  critical  Reynolds  number  baaed  on  a reduced  velocity 
within  the  boundary  layer. 

When  the  Reynolds  number,  referred  to  the  mean  ordinate  of 
the  velocity  profile,  becomes  compatible  with  stability  cf 
turbulence,  then  this  kind  of  flew  is  building  up  and  spreading 
around  sc  as  to  extend  soon  all  over  the  cross  section.  Thus, 
finally  a turbulent  agitation  becomes  superposed  over  the 
molecular  one  throughout  * 


Continuity  of  Velocity  Profile  , Exponential  Transition 

15.  In  two  different  ways,  Reichardt  (I94.G)  and  Elasr  (1914-9) 
have  expressed  analytically  tha  continuity  of  the  velocity 
profile . 

16.  In  Reference  (m^  ),  an  Intermediate  exponential  law: 

m/^  » \s.s  - tu.s)  (,  *177)  (4) 

hC 

is  interpolated  between  the  linear*  and  the  logarithmic  region 
in  a manner  giving  a continuous  transition  of  slope  at  both 
ends  of  the  intermediate  rone. 

I?.  In  Reference  (qw  ) the  continuity  of  the  velocity  profile 
is  satisfied  with  another  exponential  function: 


ij/* 

*#*•» 

\nr 

hotli  kj 

10 

0 4 *]  * **■*•  ~9 


(5) 
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18.  This  exponential  curve  is  starting  directly  from  the  wall 
(no  laminar  sublayer)  and  leading  smoothly  into  the  logarithmic 
branch  at  the  ether  end. 

19.  However*  neither  of  the  above  solutions  seems  to  satisfy. 
In  proximity  of  the  wall,  the  requirement  of  continuity  for 
the  group  of  the  fluctuating  instantaneous  components  (u,  , v#, 
w,  ) that  express  the  turbulent  motion  as  superposed  to  tne 
mean  flow. 


Relohardt ♦ s Unified  Law  for  the  Complete  Velocity  Profile 

t.  Boundary  Layer  and  Centralized  ?r- and tl -Number 

20.  More  advanced  investigations  of  recent  date  (Reicherdt, 
Reference  (m^) , 1950  > have  led  to  defining  .a  characteristic 
thickness  of  the  hsat  conducting  boundary  layer  (^J^)  as  the 
dimensionless  distance  from  the  wall  of  an  imaginary  dividing 
line.  At  this  distance  thsre  is  equality  between  the  molecular 
( k ) and  the  turbulent  ( ) coefficient  of  thermal  conduc- 

tivity. Among  the  two  causes  of  heat  conduction  at  a point 
with: 


^ : tha  molecular  one  prevails. 

•>»  > T : t^16  turbulent  one  prevails. 

21.  Similarly,  another  characteristic  thickness  (‘?}|  ) is 
defined  for  the  viscous  boundary  layer.  It  corresponds  to  an 
equality  between  molecular  ( f*  ) and  turbulent  ( A^)  coefficient 
of  viscosity.  Among  the  two  cauaea  of  shear  =tre«s  at  a point 
with: 


! the  molecular  one  prevails. 

- % 1 the  turbulent  one  prevails. 

1 "i, 

22.  Those  definition#  presume  a co-existence  of  molecular 
and  turbulent  effects.  The  preponderance  of  the  former  or  the 
latter  is  depending  on  whether  the  dimensionless  distance  ( ) 

from  the  wall  is  smaller  or  larger  than  the  characteristic  1 
thickness  of  the  respective  boundary  layer. 


2j.  Considerations  based  on  a thorough  analysis  of  measure- 
ments by  Reichardt  and  other  investigators  have  shown  that 
for  a sufficiently  large  '»’}  , the  molecular  effects  become 
negligible  as  compared  to  turbulence.  Yet  a converse  assump- 
tion of  a purely  laminar  sublayer  at  wall  proximity  is  not  free 


- f1  » 4.  *<9  4 * 4 <1 

auxv  uavuu  < 
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24,  The  layer  with  prevalent  molecular  viscosity  ( ^ ) 

-.eeas  to  have  a fairly  constant  thickness,  say  * "+  ° A* 

This  value  of  ■*),  t corresponds  toi  4.  = r.e  is  the  logarith- 
■lc  profile  (cf.  f omuls  (2),  p.  9)*  However,  the  thickness  ( ) 

of  the  heat  conducting  layer  appears  strssgly  dependent  on  the 
general  Pranfltl  numbsrji  ^ ^ ^ 


The  actual  picture  1st 
1*!  > I « — * ^ 

< l . — •»  >11, 


and  experiments  have  shown  thets 


! 4 V*  < * 

r 


25.  Ass 


. * /- 

Jl  ,'k~ 

vr  T‘/r~ 


At 


(6) 


it  follows  that  the  relative  influence  of  turbulence  on  heat 
conduction  is.  for  actual  coalitions  of  'V/4<-  > / , larger 

than  it*  relative  influence  or:  diffusion  of  aoawnttm  (viscous 
effects),  >roa  good  measurements,  Rsichardt  has  computed  quite 
recently  for  a range  of  the  general  Prandtl  numbers 
0.?2  4 4 ‘.es-o  s a corresponding  rangsi 

;i‘*  ^ 'i.  ^ for  the  thickness  ijoT  the  heat  conducting 

boundary  layer. 


26,  It  follows  that,  at  higher  values  of  Pr‘,  there  is  a large 
inter^sli 


io 
• — 


(relatlvelys  large  minuend,  small  subtrahend) 
within  which  an  assumed  laminar  sublayer  would  have  to  accom- 
modate a prevalently  turbulent  heat  flow  with  a purely  laminar 
transport  of  momentum.  This  would  be  s basic  contradiction, 
and  t WV?  native  cf  complicating  the  theory  by  implying  a 

different  mechanism  of  turbulent*  for  the  transport  of  heat 
than  ox  momentum  would  not  appear  less  objectionable.  Hence 
the  necessity  of  thinning  down  the  sublayer  so  very  orach  as  to 
having  it  dropped  completely. 


12 
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b.  Wall  Region 

27*  Reichnrdt  took  in  account  the  continuity  requirement  with 
respect  to  the  three-dimensional  group  of  the  fluctuating 
parts  (u*,  , wt)  of  the  velocity  components.  Thus  he  obtained 

(19f?l)  the  following  expression  for*  the  ratio  ( /*  ) between 
turbulent  and  molecular  shear  stress  in  the  vicinity  of  a wall t 


ft 

1 


3 

1 


(7) 


29.  In  recent  papers  (References?  (m  ) „ (blJ  ) It  is  aseumed 
that,  at  the  wall,  the  derivatives  3 


\ 

i 


are  correlated  to  each  other  so  ay  to  yield  a finite  value  for 
the  average  of  their  product* 


29«  It  follows  that  the  established  boundary  conditions  of 
non-slip  and  finite  shear  stress  at  the  wall  have  tc:  be  supple- 
mented bv  added  requirements.  Thus,  not  only  the  function 
( Arif*’  ; itself,  but  Its  first  and  second  derivative  as  well, 
with  i-ogord  to  -'i  , hove  to  vanish  zt  the  wall. 


30.  Rigorously,  such  a prevalence  of  the  cubio  term  of  a power 
series ; 


is  well  established  only  for  a partly  developed,  two-dimensional 
state  of  turbulence  close  to  the  wall.  Yet  this  is  to  hold, 
most  likely,  as  well  for  the  state  of  the  fully  developed, 
three-dimensional  turbulence  (Reference  (m^  )). 


4 

U 


f 


1 
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31.  The  integrated  set  of  four  boundary  conditions  at  the 
wall  and  an  assumed  asymptotic  trend  of  the  function  F>  toward 
the  linear  form  vCy  at  some  distance  from  the  wall  constitute 
a solid  physical  background  for  the  new  Reichardt  function 
(Figure  2) 1 


* w C >)  - 1.  'h, ) 

32.  For  a flow  with  drop  of  pressure,  the  identity t 


(8) 


may  bo  easily  transformed  into: 


a, 


\ * 


(?) 


(10) 


where  the  expression  '/t+p  approximates  fairly  well  a flow 
close  to  the  wall  of  prismatic  ducts*  The  croaa  Faction  may 
be  either  closed  (tube)  or  open  (channel).  This  was  obtained 
with  a linear  drop  of  shear  stress  from  its  value  ( 'c'„  ) at  the 
wall,  and  the  neglect  there  of  the  term  ( y/n.)  as  compared  to 
unity. 


33.  The  insertion  of  function  (8)  into  (10)  lends  to  the 
following  differential  equation: 


a* 


( t + k>i  \ f \ — 


" ii 


u Vy) 


*5  Ii 


\ 4K*) 


irxuojraoiOu  u y pares  yxoxas: 


(10a) 


♦ 


-i-  Xw_  f > -v 
K 


+ 


K~nt  

w O ’ • *1 


(11) 


35*  Unfortunately,  the  integral  ( ^ ) of  the  right-hand  side 
turns  out  to  be  intractable.  Approximate  solutions  ha^e 
developed  (Reference  (m  ))  in  two  ways;  these  are: 


1* 


* 1 
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a,  G/raphical  integration  end  negl&ot  of  the  unity 
against  K}  for  higa  values  of  -tj 

b.  Replacement  of  the  intractable  integrand  by  an 
approximate  substitute  such  as  to  make  the  integration  possible* 
Prom: 

1 ^ * a U,  h<, -.)wt-d)|  * £ 

1 i di 


It  follows  that : 


£ s C,  1 I - -5-  l-mj 

* | * l 1 * ' 1,  • J 


36.  After  carrying  out  the  Integration  and  determining  the 
relevant  constants  from  the  boundary  conditions,  Reich&rdt 
obtained  finally,  for  the  wall  region,  the  velocity  profile: 


'P  =.  “-ii*  ~ * c,  j \ (12) 

with  the  constants  ( cf.  p.  9 for  5s 

K x ‘/C,  * ( 

c (-  i* 

i - o < 

& CS  */» 

37.  In  Figure  3 the  solution  naw  (graphical  integration)  is 
marked  by  a full  line  whereas  the  dashed  one  represents 
solution  "b"  (approximate  integration).  Between  "as  and  "b”, 
only  lnsigyiiricant  discrepancies  with  alternating  sign  occur 
for  certain  regions  of  the  dimensionless  distance  ( *t  ) from 
the  wall. 

In  bha  same  figure,  the  ratio  x±(v  is  plotted  versus  ; 

> 4.  i-  . M _ i JO 1 i.*  * 


1 - M _ A JO -1  *-  ^ 1 AA 

uaao  xaubvi'  as  ran^niK  a aoaii  a ou  av^u 


c.  Median  Region 


38s  In  the  median  region  of  a flow  through  prismatic  ducts, 
the  agitation  is  prevalently  turbu?.ent  I.e.  the  effects  of  the 
molecular  motion  on  ths  velocity  profile  may  be  neglected. 

The  distance  ( ^ -=*  n-'f  ) from  the  median  plane  (channel),  or 
the  axis  (tube),  is  a more  convenient  coordinate  than  y. 


~ai- 


T 

i 
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Therefore,  the  orgin  of  the  ordinate  in  a direction  normal  to 
the  wall  is  shifted  from  the  wall  to  the  center. 

39*  In  Reference  (m^),  the  relevant  velocity  profile  is  based 
on  accurate  empirical  data.  The  relative  coefficient 


of  ths  apparent  shear  stress  was  obtained  by  directly  measuring 
the  derivative  of  the  stagnation  pressure  throughout  the  cross 
section  with  a double  Titot  tuba.  Such  experimental  values 
obtained  in  a channel  of  rectangular  cross  section  (Figure®  ii 
and  5)  were  expressed  by  the  formula: 


(13) 


lj-O.  This  relation  is  inserted  in  the  counterpart  of  equation 
O)  which  for  the  median  region  is; 


av  . - ^ A ( ^ ) 

4 -1.  After  integrating  and  determining  the  relevant  constant, 
the  velocity  profile  becomes  for  nearly  the  entire  cross  section 
(except  in  vicinity  to  a wall)  : 


UL  -K 

-v. 


where : 


- K \r^sl 


lp2.  In  the  graph  (Figure  ij.)  the  measured  relative  velocity 
? ss  **■}**> as  well  as  its  derivative  tfiii  , and  the 
computed  second  derivative  A*y/4a‘'  are  all  plotted  versus  the 
distance  (z)  from  the  median  plane  ( z.  » y ). 


s l 
1 1 


I 


NAVORB  Report  2958 

I4.3.  It  J s Interesting  to  nota  that  9 maximum, 

halfway  between  wall  and  median  plane,  at  a station  where  the 
curvature  presents  a minimum  (Figures  I4.  and  5)  • 

Vf.  With  regard  to  a logarithmic  profile  of  velocity,  the 
measurements  were  consistent  not  only  in  vicinity  of  the  wall 
but  further-  away,  toward  the  mad  I dii  region  of  the  cro  ss  section. 
Nonetheless,  the  clot  ■-*/**  versus  log  '*7  showed,  for  still 
higher  values  of  ^ , a distinct  swerving  away  of  the 

empirical  values  from  the  straight  lino  (Figure  6).  Reichardt 
succeeded  (Reference  (m*  ))  to  eliminate  such  discrepancies  by 
showing  that,  instead  of  x/**  , a modified  expression  *-¥t /■**? 
is  a function  of  log  ■*»  . The  additional  velocity  x.,.  for 

the  inner  region  of  flow  (hence  the  subscript  1)  was  defined  by 
the  function  (Figure  6} 1 


« vf  ’ 


uf(»»VJ  1 
» 4 tlV.^1 


d.  Complete  Profile  All  Over  the  Cross  Section 

!^5*  If  the  expression  (15)  for  is  entered  in  the  universal 
logarithmic  lawi 


^ , * 


then  the  following  unified  law  is  holding  in  esse  of  turbulent 
flow  for  the  velocity  profile  of  the  entire  cross  section: 


The  complete  velocity  profile  of  turbulent  flow  in  • cylindrical 
tubs,  is  shown  in  Figure  7 for  different  Reynolds  numbers. 
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4.6.  It  ia  easy  to  verify  that  the  general  relation  (17)  con- 
tains as  particular  cases  both  the  velocity  profiles  that  were 
previously  derived  for  the  wall  region  and  the  median  region 
respectively.  The  former  (equation  (12))  occurs  for  ><. i ,x  « •' 
The  latter,  corresponding  to  large  values  of  >j  . takes  the 
form : 


aa 


x.r(  \ t V<)  I 

— 1 f + < 


(18) 


Itf . The  constant  term  ci  is  expected  to  depend  on  the 
shape  cf  the  cross  section.  It  should  be  determined  by  experi- 
ment for  the  existing  conditions  of  turbulent  flow. 


!|8.  The  unified  velocity  profile  (17)  of  the  prismatic  duct 
is  expected  to  hold  as  well  for  the  turbulent  flow  along  a 
flat  plate.  The  latter  agrees  well  with  the  former  in  proximity 
of  the  wall  where  the  distribution  of  velocity  depends  mainly 
on  the  shear  stress  velocity* 


'/* 


I4.9,  In  regions  afar  from  the  wall,  however,  individual 
features  of  the  velocity  profile  may  appear.  The  flow  along  a 
plate  requires  higher  additional  velocities  ( 44’  ) at  the  outer 
region  than  those  used  at  the  Inner  region  of  a duct  for  a 
truly  logarithmic  profile  satisfying  the  complete  set  of 
boundary  conditions. 


18 
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Part  II : The  Profile  of  Mean  Temperature 


Introduction 


50.  This  la  the  second  part  of  a review  on  the  characteristics 
of  the  turbulent  flow  with  heat  transfer  in  a prismatic  duct  of 
symmetric  cross  section  or  in  the  boundary  layer  of  a flat 
plate. 


51.  For  th6  laminar  flow  with  such  a geometry,  the  temperature 
profile  has  been  rigorously  computed  some  wb^le  ago.  Yet,  for 
the  thermal  field  of  turbulence  r.o  general  analytic  solution 

is  known. 

52.  A main  objective  of  this  part  is  the  recent  progress  in 
solving  the  turbulent  case  by  integrating  the  temperature 
equation,  numerically  or  graphically,  in  successive  steps. 

53.  The  turbulent  flow  is  taken  as  fully  developed,  and  aB 
feeding  upon  the  energy  of  a forced  mean  motion  that  is  assumed, 
stationary. 

5I4-.  The  solid  boundary  ia  postulated  aa  smooth  and  of  nearly 
uniform  temperature. 

55*  It-  is  found  by  integration  how  the  fluid's  temperature  is 
dealt  out  normal  to  the  wall  in  terms  of  either  the  distance 
from  it  or  the  velocity.  The  latter  argument  seems  preferable 
for  the  higher  Prandtl-numbers . 

56.  The  variables  are  made  dimensionless  on  a number  of 
different  bases.  The  divisor  is  either  one  of  the  border 
values  (at  wall  or  free  stream)  of  the  relevant  quantity  or  its 
mean  value  ova  I*  tliW  boundary  layer. 

5?.  The  molecular  and  turbulent  part  of  both  shear  stress  and 
heat  flux  density*  are  expressed  in  a function  of  the  gradients 
of  velocity  and  temperature  respectively.  Such  combined 
expressions  may  lead  to  virtual  solutions  for  the  temperature 
in  the  integral  form. 

58.  The  integrands  depend  upon  the  fluid's  influential 
properties:  f and  assumed  for  the  time  being  as  nearly 

constant  - and  on  three  func^onal  arguments.  As  such,  the 
following  dimensionless  expressions  are  conveniently  chosen 
(References  (in)  and  (q))» 

* flux  = time  rate  of  flow;  the  flux  density  is  taken  per  unit 
of  cro^s  area. 

19 
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a*  The  ratio  * /*  » *7  /*  between  the  turbulent  u*rt  of 

tha  ahaar  ■ trass  (flux  density  of  mo  man  turn)  and  its  oulwouiar 
counterpart.  T'ha  exact  form  ot  * A B haa  baan  discovered 
recently  by  tha  inductive  method  (References  (at_).  (ay), 

I 1 P * 

im« ) , • 

* 

b.  Tha  dimanaionlaes  heat  flux  density*  ^ , based 
on  Its  value  at  the  vail*  The  function  may  bo  calculated 

in  tarns  of  the  dimensionless  wall  distance  (Reference  (■-) } by 
integrating  the  equation  for  the  continuity  of  the  heau  flow, 
provided  the  shear  distribution  la  known* 

e.  A function  (&)  to  account  for  the  experimental  faot  - 
inf erred  from  References  (g)  and  (i)  - that  turbulenoa  is 
diffusing  heat  faster  than  momentum.  This  la  expressed  by  the 
ratio  >rv  • A%.  j /v  between  the  coefficients  ; " A*  " for 
the  thermodynamic  effect  of  turbulenoa  and  ” A..  " for  its 
mechanical  action  as  apparent  Yiacoslty. 

59-  Tha  form  of  this  important  function  ia  unknown,  though 
ita  boundary  values  are  pretty  sure  ( I £ -m.  < Z ),  and 
average  values  over  the  boundary  Isjers  nave  bean  estimated 
(Figure  1)  for  particular  cases  (prismatic  duct,  plate). 

60.  Under  the  above  premises,  the  integration  for  the  profile 
of  mean  temperature  haa  been  recently  performed  in  successive 
steps*.,  There  la  c rapid  convergence,  especially  for  not  tco 

<M«e1  1 Osa_«Miw\t«iAA 
■SSSSAA  a i -UUWKV4  Q t 

61.  By  numerical  or  graphical  integration,  a satisfactory 

approximation  la  reached  in  one  or  two  steps  for  air  and  fluids 

with  higher  Pr-number n , while  three  or.  at  most,  four  steps  are 

needed  in  the  case  of  very  small  ?r -number a (References  im  ) . 

# % • \ * ' 

62.  The  temperature  field  is  discussed  in  view  of  the  recently 
extended  e«i «ui scions  of  the  heat-  transfer  over  a range  of: 

e.  Very,  high  values  of  the  Pr-number  (References  (f)  and 
(m) , up  to  t -•  ) • 

b.  Very  low  Pr-numbers  (S«xerono«i  (a),  (q)  end  (s)  down 
**°  T*  ■*  mT*,  I®”*  I O- *,  ) • 

63.  The  former  category  corresponds  to  son-metallic  liquids 
like  water  and  oils.  The  latter  covers  conditions  of  motele 
in  liquid  state  (mercury,  molten  metals  and  alloys).  The 


* Tapers  with  analytical  solutions  for  the  single  particular 
oaa <m  Pr*l,  m = l,  and  independency  of  the  fluid’s  proper^ 
ties  from  temperature,  are  not  considered. 
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lower  range  become  particularly  interesting  of  late,  in  view 
of  the  necessity  of  removing  heat  in  large  amounts  from 
nuclear  furnaces,  at  a high  temperature  level,  by  applying 
molten  metals  as  coolants. 

64.  The  shape  of  the  temperature  profile  appears  to  be  highly 

dependent  upon  the  generalized  Pr-numbsr:  f<J  . ?%. 

(Figures  9,  10  and  11),  whereas  it  seems  less  affected  by  the 
Reynolds -number . However,  the  dimensionless  heat  flux  density 
i»  found  more  dependent  upon  the  Re-numbor. 

65.  Almost  irrespective  of  both  parameters  (Pr'  and  Re),  the 
average  mean  temperature  for  turbulent  flow  in  a circular  tube 
is  found  at  a distance  equal  to  0.3  of  the  tube  radius  from 
the  wall  (References  (p  ),  (p  ),  Jakob,  M.  and  (q  ),  Elssr,  K„) 

*»  'A 

66.  By  a sustained  difference  of  temperature,  between  a solid 
wall  and  a fluid,  a heat  flow  is  caused  between  them  in 
accordance  with  the  formation  of  a thermal  boundary  layer = The 
steady  state  of  such  a field  of  enthalpy  and  temperature  is 
investigated  for  generally  specified  geometric  conditions  at 
the  boundaries. 

67.  In  addition  to  the  external  cause  of  heat  flow,  there  is 
the  viscosity  from  molecular  and  turbulent  origin  to  account 
directly  or  indirectly  for  internally  released  heat  from 

A i <*  !d1  AO  til  rm  r>*f  rnA^Vianl  nal  UnArcnr  vrV**ri  fl  t « 

sufficiently  rais  6 Ci . 

68.  The  thermal  field  is  governed  by  the  fluid's  Pr-number, 
in  addition  to  at  least  the  Re-number-  of  the  flow,  as  long  as 
the  compression  is  not  so  high  as  to  involve  the  Mach-number  as 
well.  Anyhow,  the  compressibility  effects  are  presumed  to  be 
almost  negligible,  Moreover,  neither  combustion  nor  changes 

of  phase  are  postulated.  The  latter  restrictions  allow  for  a 
similarity  between  the  fields  of  enthalpy  »;•.;]  t- sup s r s. tur 0 . 

69.  Since  the  early  work  of  Osborne  Reynolds  (I87I4.,  189c), 
the  effects  of  turbulent  agitation  on  diffusion  of  momentum 
and  heat  are  currently  expressed  by  analogous  (Figure  1) 
formulas  (G.  I.  Taylor ; Prandtl.  von  Harman)  to  those  of 
corresponding  molecular  effects  (Newton,  Fourier).  The  flux 
density  of  both  momentum  (shear  stress  x*  ) and  heat  ( % ) are 
split  into  a molecular  part  and  a turbulent  counterpart. 


X 
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70.  Th=  respective  contributions  are  expressed*  in  function 
of  coefficients  and  cf  gradients  of  velocity  end  temperature 


as  follows: 

A 

- -rv* 

r 

t 

(19) 

1 

f 

I! 

1 

- crf  Yt 

71.  At  the  wall  (subscript  0)  the  turbulent  agitation  is 
completely  impeded.  This  means  A a.o  * A • Hence* 


(19a) 


i£m 


72.  Above  expressions  (19)  are  conveniently  combined  in  a 
double  ratio : 


d i- 

r*tK 


it  it  .vf, 

at  a 


ir 


_ 'll  Pr 

V(  iT 


(20) 


73.  By  the  dimensionless  quantity  "Pr’*  or  generalised  Pr- 
number,  the  fluid’s  thermal  property  "Fr"  is  linked  with  a 
characteristic  ratio  "m"  between  two  effects  of  turbulence. 
Therefore  RPr’"  is  2.  governing  par  excellence  ror 

the  turbulent  flow  when  transport  of  both  energy  (heat)  and 
momentum  occurs. 

7^-.  Through  turbulence,  heat  seems  to  spread  faster  than 
momentum.  It  is  a fact  concluded  from  extensive  experimental 
evidence  (References (d) , (g),  (m) , (s)). 


* henceforth  unFarred  symbols,  i,  T,  u,  v,  etc.  are  used  when- 
ever they,  unmistekenly , refer  to  mean  temporal  values,  unless 
it  is  otherwise  expressly  stated.  However,  by  necessity,  the 
vincula_ax*6  maintained  for  denoting  an  average  of  a product. 


say 


f * 


Wt 


T*  rjr 

f i 
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75.  a measure  of  this  phoneme non  is  the  ratio  "m"  betwasn  the 
coefficients  "A."  for  tha  thermodynamic  effect  of  turbulence, 
and  n for  i£s  mechanical  action  by  apparent  viscosity.* 

?6.  The  ratio  of  these  coefficients  - ss.  V/V  is  a 

characteristic  function  of  the  share  taken  by  ths  turbulent 
effects  with  regard  to  the  heat  flux  density  and  to  the  internal 
release  of  heat  by  viscosity.  As  such  "mn  seem*  to  be  dependent 
on  the  wall  distance,  and  Pr*  varies  also  with  the  local  Prandtl- 
number * The  latter  changes  with  temperature  - except  fer  ideal 
gases  - and  the  dependency  is  particularly  strong  for  viscous 
fluids. 

77.  By  neglect  of  the  thermal  variation  of  wPr",  the  study 
would  be  limited  to  small  changes  of  temperature  for  these 
fluids  whose  properties  strongly  depend  upon  it*  Yet,  by  using 
an  averaged  Pr  (over  the  boundary  laysr) . some  improvement 
would  seem  attainable. 

78.  ‘ No  generally  accepted  theory  about  the  mechanism  of 

turbulence  seems  yet  available  for  yielding  the  form  of  the 
function  '’m"  over  the  viscous  and  thermal  boundary  films.  Yet, 
it  is  pretty  sure  that : I < y*.  < 1 j viz.  "xnn  has  the  border  value 
(1)  at  the  wall  (0)»  -rvio  « \ , while  at  the  free  stream  end 

of  the  boundary  iayer  it“ seems  to  approach,  » without  quite 
reaching  - the  value:  -m.  *>  *■  of  the  free  turbulence  (hot 
Jet  or  wake  of  heated  rod) . 

79.  As  the  shape  of  the  function  “ra”  is  unknown,  ths  estimated 
average  values  (Figure  l)  over  the  boundary  films  are  rather 
uncertain. 

80.  The  measurements  of  temperature  and  velocity  profiles  for 
air,  flowing  over  heated  bodies,  seem  to  show  a dependency 
between  the  average  value  of  nnr  over  the  boundary  layer  and 
the  layer's  geometrv.  It  matters  whether  the  flow  occurs  in  a 
duct  ( channel f tube)  or  around  a submerged  body  (flat  plate, 
missile,  wing;. 

81.  It  may  be  of  interest  to  examine  whether  - and  to  what 
extent  • flow  conditions  such  as  the  wall  roughness  ard  the 
pressure  gradient  (positive,  negative  or  null)  influence  the 
distribution  of  ”ro" . 

82.  In  order  to  find  the  heat  transfer  at  the  wall,  the 
velocity  profile  is  needed  In  general.  ‘Similarity  may  exist 


«•  Both  A*  and  , called  also  "Austausch"  coefficients, 

have  the  dimensions  of  the  dynamic  coefficient  ( j*  ) of  viscosi  ty  <. 
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betw&en  the  profiles  of  velocity,  shear  stress,  and  temperature, 
regardless  of  the  Re-number,  insofar  as  the  finite  boundary 
layer  Is  laminar  or  the  flow  occurs  at  a Mach  number  high 
enough  to  reduce  the  width  of  the  boundary  layer  to  the  order 
of  a molecular  mean  free  paths 

83.  On  the  contrary,  for  turbulent  flow  in  subsonic  or  low 
supersonic  range,  the  thickness  of  the  viscous  boundary  layer 
is  strongly  dependent  on  the  Re-number. 

8h..  No  yielding  of  the  temperature  profile  may  ensue  by 
necessity  from  knowledge  about  the  velocity  profile.  In  general 
a separate  integration  should  be  indispensable  for  obtaining 
the  former  near  the  wall.  This  point  would  appear  to  vindicate 
devoting  to  the  wall  region  of  the  thermal  field  a main  part 
of  our  survey. 


85.  Moreover,  it  is  important  to  know  the  temperature  profile, 
because  the  heat  transfer  at  the  wall  depends  substantially  on 
its  shape.  The  heat  transfer  coefficient  at  the  wall  1st 


vMS) ml 


86.  For  smaller  Pr-numbers  the  temperature  function  seems 

more  convenient,  while  for  higher  Pr-numbers,  the  function  {m.) 
appears  preferable.  It  Is  interesting  to  note  that,  over  & * 
limited  range  of  the  Pe'clet  number  ( Re  Pr  > 2500),  the  heat 

transfer  might  as  wall  be  derived  from  a general  relation  not 
linked  to  knowing  expressly  how  the  temperature  is  dealt  out 
(Reference  (mr)).  The  letter  possibility  will  be  discussed  in 
the  third  part  of  the  review  dealing  specifically  with  the  heat 
transfer  at  the  wall. 

87.  Insofar  as  the  Influential  properties  of  the  fluid  (c„.k. 

■ft  and  o ) are  nearly  constant,#  the  energy  equation  wilV 

remain  linear.  Therefore  a superposition  should  be  permissible 
in  regard  to  the  thermal  fields  of  external  and  Internal  origin. 
The  latter-  is  baaed  on  adiabatic  walls  (perfect  thermal 
Insulation)  at  equilibrium  temperature. 


* A dependency  upon  temperature  - although  weak  - would  seem 
more  justified  for  the  specific  heat  than  fer  the  heat 
conductivity. 
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83.  With  turbulence  included  in  the  aero-thermodynamic 
problem,  there  is  an  involved  degree  of  intricacy  in  the 
combined  flux  cf  momentum  and  heat. 

89.  Three  kinds  of  motion,  basically  different  in  Huturd  | Mi  G 

interacting,  i.e.,  the  forced  mean  motion  and  two  sorts  of 
agitation  on  different  levels ^ One  of  them  operates  on  the 
tiny  mcleoular  scale  by  transport  of  momentum  (viscosity)  and 
heat  (conductivity),  while  the  turbulence  occurs  on  a molar 
scale.  Its  relevant  effects  are:  apparent  viscosity  and  heat 

convection. 

90.  Under  such  involved  conditions,  no  analytical  integration 
is  in  sight  for  the  partial  differential  equations  that  govern 
the  conservation  of  momentum,  mass  and  energy  in  the  general 
turbulent  flow  with  heat  transfer.  Neither  are  general  solu- 
tions attainable  yet  for  the  simpler  case  of  a plane  mean 
motion  (the  turbulent  one  is  always  three-dimensional),  as  that 
occurring  afar  from  the  ends  of  prismatic  ducts  or  along  a 
flat  plats. 

91.  However,  with  simplifications  that  wars  moderate  and 
Judicious,  tho  thermal  field  was  integrated  numerically  and 
graphically  for  the  turbulent  flow  in  a straight  circular 
tube  (References  (m)  and  (q)). 

92.  In  regard  to  a flat  piste,  a similar  procedure  seems 
promising  only  at  a very  large  Re-number,  for  it  allows  the 
reduction  of  the  plane  temperature  field  to  a mere  one- 
dimensional problem.  Other  attempts  to  convert  the  two- 
dimensional,  partial  differential  equation  into  an  integrable, 
ordinary  differential  equation  were  shown  to  be  not  rigorous 
enough  (Reference  (q  )}r 


Underlying  Assurrptions 

93*  To  recapitulate,  the  present  study  is  postulating: 

a.  A fully  developed  turbulenc  motion  feeding  upon  the 
energy  of  a stationery,  forced  mean  flow  in  a prismatic  duct 
or  along  a fist  plate. 

b.  A smooth  solid  wall  throughout. 

c.  A heat  flow  caused  by  a difference  of  temperature 
between  the  solid  wall  and  the  fluid  in  the  free  stream  in 
accordance  with  the  3hape  of  the  thermal  boundary  layer. 

d.  A practically  constant  wall  temperature  over  the 
entire  surface  of  the  solid  boundary.  Accordingly,  tho  trans- 
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port  of  both  heat  and  momentum  in  the  respective  boundary 

Ai  <S  ..  . » ..  **  *-  . M.  « 1 

1 iiiua  J_a  HUA'i'itlx  VU  OlAO  noxx. 


e.  Flow  conditions  such  as  to  not  affect  appreciably 
the  density  and  other*  influential  properties  of  the  fluid 

< . f-  . >• 

f.  Negligible  effects  of  radiation  on  the  heat  economy.. 

94*  By  postulate  d,  the  resultant  heat  flow  becomes  one- 
dimensional in  a convenient  direction  for  making  tractable  the 
differential  equations. 

95.  It  would  seem  advisable  to  promote  clarity  and  generality 
by  not  trying  to  write  out  too  soon,  under  the  integral,  the 
imperfectly  known  functional  expressions  for  /*•**■/?*■  and 

•wv.  * Ay  / • Besides,  as  to  the  function  "m" , as  well  aa 

to  the  influential  properties  of  the  fluid  (P r, />  ).  only 
average  values  may  be  anticipated  adequately,  for  the  time  being. 
Later,  the  variable  functional  arguments  of  the  mixed  flow  and 
of  the  fluid's  properties  might  be  taken  into  account,  provided 
they  are  reducible  to  a Judicious  form  and  an  acceptable  one  for 
keeping  the  integrals  tractable  without  oversimplifying.  It 
seems  important  to  try  to  reach  nontrivial  solutions  by 
non-too-arduous  method a. 

96.  The  above  set  of  postulates  proved  apt  to  render  soluble 
the  differential  equations  for  the  turbulent  flow  with  heat 
transfer  by  numerical  or  graphical  methods  of  successive 
approximation  (References  (m^),  (m r) , (q^)).  The  convergence 
proved  rapid  enough  to  yield* satisfactory  answers. 


Relchardt ' a Refinements  about  the  Turbulent 
Boundary  Layer 

97.  A main  incentive  to  the  recent  developments  In  the  boundary 
layer  theory  was  the  trend  toward  approximating  closer  the 
temperature  profile  for  turbulent  flow  and  extending  the 
solutions  over  a much  wider  range  of  Pr-numbers.  In  this 
connection  much  attention  has  been  given  of  late  to: 

a.  The  boundary  conditions  for  both  the  mean  motion  and 
the  superposed  turbulent  one. 

b.  The  observed  continuity  of  decay  of  turbulence  from 
the  free  stream  to  the  Immediate  vicinity  of  the  wall. 

98.  The  early  partition  of  the  boundary  layen#  in  on©  purely 
laminar  part,  adjacent  to  the  wall,  and  one  purely  turbulent 
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outer  part,  allowed  or.ly  n single  value  of  the  Pr -number ; 
viz.  Pr  ~ 1 with* 


0 , in  the  laminar  sublayer 

/ 

/ 

^*1  , in  the  turbulent  layer 


99.  Besides  this  limitation,  such  a dichotomy  of  the  boundary 
layer  implies  a severe  discrepancy  in  the  form  of  a noncbservable 
singularity  of  the  temperature  profile  at  the  outer  limit  of  the 
purely  laminar  sublayer. 

100.  Later  - the  attempt  to  deal  with  heat  transfer,  by  means  of 
fluids  with  Pr  / 1,  led  to  the  well-known  and  still  rather 
popular  trichotomy  of  the  boundary  layer.  This  meant  inserting 
a buffer  layer  (simple  or  composite)  of  a mixed  - laminar  end 
turbulent  - nature,  between  the  laminar  sublayer  end  the 
turbulent  outer  reg  Ion m Hovsvcr  km  - A_  =1  was  maintained 
in  the  two  outer  subdivisions,  and  "the  tripartite  boundary 
layer  did  not  escape  from  ell  the  previous  shortcomings  by 
failing  to  satisfy  certain  boundary  conditions. 

101.  To  get  rid  of  such  contradictions,  is  the  aim  of  the 
newest  phan*  in  the  development  of  the  boundary  layer  theory. 

It  means  abolishing  the  artifice  of  a purely  laminar  sublayer 
and  inferring  a refined  distinction  between  a characteristic 
width  ( ^ ) for  the  thermal  profile  and  a characteristic  one 
( ~'j,  ) for  the  viscous  boundary  layer.  Except  for  pr  : * 1, 
the  boundary  layer  shows  a different  characteristic  width  in 
the  fields  of  temperature  and  velocity. 

102.  The  width  of  a film,  in  which  a molecular  effect  is 
prevalent  over  Its  turbulent  counterpart,  is  expressed  as  a 
limit  value  of  a dimensionless  distance  from  the  wall,  (say  >/t 
or  preferably*  *»>  «.  ?*/*«•) . Within  a so  defined  edge  of  the 
velocity  field,  Newtonian  viscosity  prevails  over  the  apparent 
one  due  to  turbulence. 

Similarly,  for  the  temperature  field,  a character is tic  thermal 
layer  at  the  wall  marks  a region  of  prevalence  for  the  Fourier 
conduction  ox  heat  by  taolecu^ar  agitation  against  the  thermal 
convection  from  turbulence. 

103.  With 
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the  generalized  Fr-rtmber  Is:  «*  *h  % • '//%  . Among  tha 

four  functions:  a,  P r,  / and  , only  the  last  mentioned 

( Ar/t*-  ) appears  a*  sufficiently  explored  ©Tar  the  boundary 
region  {Rsfsrasicea  (m)  and  (q)},  Howerar,  with  rough  estimates 
for  "s"  and  "Pr",  the  function  "a  ^ may  be  approximated  froan 


j - } {y* ) /*{*!} 

'm  f 

loif.  Th«>  sxact  form  of  " A * as  a continuous  function  of  * *}  " 
as  discovered  recently  (References  (mA) , (mg) . (a/)}  by  tbs 
inductive  method  is  baaed  on  tha  bast  available  measurements 
and  real  boundary  condition*.  It  also  aatisfiae  more  oompletely 
tha  continuity  requirements  for  both,  tha  mean  motion  and  tha 
turbulence. 


105.  With  the  d^sionl...  argument  wither  in  the  form  of  a 
Ha-number  (•»**.  ) of  (j  ■>  1/*)  for  tha  wall  region  or  in 

the  form  %/«.  for  the  core  ( tube 'or  duot)  and  the  constants: 

K " 0.1;  and  lb  & £MA&h*  following  lsr=  appear  to  satisfy  tha 
foregoing  requirement* / 


106.  A -Region  with  **}  < 6 

Three  term*  arc  contributing  to  tha  function  fi  ae  follows: 


A 1a.  A 

rr>  1 


- » <"M 


) c.  c "t* u*4  ff 


'K  . 


1*0 


Variables : j*  t 


However,  for  extremely  small  Talus a of  „ th*  power  series 
would  saass  to  start  with  a term  proportional  to  *j’r  . This 
means:  C-s  £ 0 t cf<»  0 . Than  the  abtrva  law  would 

approach,  according  to  Referenoa  (&  )t  for  Ttry  small  , 

tne  fora: 


-r  _r 

A et  l*\  15  ■*! 

i • 


(22a) 
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107.  B-Region  with  * < ^ < »a  . 

Two  torms  only  would  contribut#  to  tha  function  f*  ( "* 7 ) 
S5  follows  ? 


5 ( c,  r V 1 

i»o,  1 


Kh  - ^*.***4%  ] (23) 

v. 


Variables  t t S'  . 


1 

1 

HT 

0 

l 

^ 1 

I 

\ 

\ *‘/ir 

"Tl 

■ 

Constants: 

K.  *: 

O.V 

II 

109* 

C-Region 

with 

■>1  : 

• 

> 10 

r 

Tills 

la  the  core  of  a 

tube 

or  the 

median  region  of  a channel. 

Bore,  " /*  " depends  not  only  upon  tha  dimensionless  distance 
( 'h/'i  ) from  the  tuba-axis  or  the  channel* • medium  plane, 

but  also  on  the  Re-number,  Indirectly,  through  the  quantity* 


*». 


f 

* K« 


7<i)  \1 


«AM. 


109*  Tha  form  of  the  function 
(m< ) , ( t) ) ia  here. 


m «?i 

..  , 


(References  (m  ) , (mv). 

^ • 


,*  H 


* Ki*  i » ^ 1 7c.\  i i t - vV**)  1 

r 5“  L * " i i I 


rp).  ^ 


110c  with  regard  trs  th§  effects  of  viscosity,  s chsracte'’'^  t(n 
layer-wiuth,  " 7,  M measured  from  the  wall,  may  be  defined  for 
the  prevalence  in  it  of  the  molecular  over  the  turbulent  effeet- 
{ 71  ^ > Tr  )•  Thus  ” 7J,  ” ia  that  value  of  the  argument 


for  which* 


f>c7  J 


(25) 


Ills  Similarly,  another  characteristic  width  ” 1*"  may  be 
referred  to  the  heat  flux  density  in  the  boundary  layer.  It 
marks  a limit  distance  from  the  wall  for  prevalence  in  it  of 
the  molecular  conduction  over  the  turbulent  convection  of  heat 


«-7 


that  value 


KAVORD  Ecpcri  2958 

( ^ >>  ).  ?hui  " *»]  " is,  for  a I, 

of  >j  for  which  1 '** 


I 1,4  ’ ?**v 


/ 5 


6) 


112*  If  Re  ia  high  enough,  then  H ends  within  the  outer 
part  of  the  B~reglcn  (large  , hence  tahh  ■*»/,,  ^1)*  Now 

rj  is  daterrsined  from  the  equation*  *» 


Vicvs-s.  : 


(26a) 

I 0«X,«.  fw 

^ t ^ "R*  | 


113s  It  follows  that,  in  general,  by  virtue  of  (25)  **id  (26): 


*■'«  < ' — \ ‘ t 

^ tut  *■  ^ 

*c  >*  — * V 

lilt*  With  7^.*  7.15  for  the  extrapolated  A-region,  the  value 
'j  * T,  * 1O.7  satisfies  the  limit  width:  l.  Very 

nearly  the  seme  characteristic  width  for  the  viscous  layer  — 
viz.  ‘*i,  » 10,8  — > is  ettaiiicd  by  using  the  ,1  -law  of  th» 
B-region  with  **|  *11* 

115*  For  a given  shear  stress  distribution,  V is  nearly 
constant.  On  the  contrary  the  characteristic  width  « ^ ) 

of  the  thermal  layer  is  variable  as  dependent  on  Pr*  ■ ^ Pt  * 

116.  It  is  interesting  to  note  that  for  Fr!  1,  the 
corresponding  low  values  of  the  function  f*  are  obtained  with 
accordingly  small  values  of  ^ *»),  . Tnat  means  an 

extremely  thin  characteristic  width  for  the  heat  conducting 
layer  with  a prevalence  of  turbulent  heat  convection  within  the 
layer  ( *1  iO  of  prevailing  molecular  viscosity* 

11?.  Example  1 for  Pr*  “ $00  and  ft  «=  0.01:  3 * P Pr*  =5, 

*•••  f although  Tt*  (0.01)  . 

118.  This  illustrates  clearly  the  fallacy  of  conceiving  a 
purely  laminar  sublayer  of  finite  thickness. 
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119.  On  the  other  hand,  for  ?r*  1,  the  corresponding 

high  values  of  one  function  /%  lead  to  7*  >j;  .Thus  the 

heat  transfer  by  metallic  fluids  operates  mainly  by  molecular 
heat  conduction  even  at  values  of  at  which  the  fjux  density 
of  momentum  is  mainly  a turbulent  phenomenon* 

120.  Example:  for  Pr'  = -'..01  and  = 1 : £ «=  (* Pr»  * 0.01 

i.e.  9. . ■ {100)2*  at  the  limit  (A^  « h ) of  the  viscous 

layer.  This  would  hold  for  a me  tallic ‘liquid  like  mercury  as 
heat  transferring  medium, 

121*  The  early  assumption  of  a common  mechanism  of  turbulence 
to  account  for  both  thermal  ( ) and  mechanical  ( ) effects 

implied  as  well  ( w ) * 1.  Ey  such  a theory  the  dichotomy 

of  the  boundary  layer  in ’ a purely  laminar  and  a purely  turbulent 
outer  part,  appeared  satisfactory  for  a fluid  with  Pr  * 1.  That 
would  be  steam  or  a polyatomic  gas,  at  least  quadratomic. 

122c  However,  according  to  the  present  views:  *• V , for 

the  flow  near  the  wall  of  a prismatic  duct.  It  follows  that 
the  restricted  validity  of  the  old  theory  deserves  being 
interpreted  in  the  light  of  the  latest  CCi  tCOp  tions.  These  show 
the  old  theory  holding  not  for  a fluid  with  Pr  = 1,  but  for 
one  with  Pr'  = 19  implying:  Pr  «=  '/l.y  » 0.72.  This 

fluid  is,  of  course,  air  and  not  steam. 

123.  The  solid  background  of  the  new  theory  is  the  continuity 
in  the  decay  of  the  turbulent  effects  from  the  free  stream  to 
the  wall,  with,  in  general,  I,  ^ Ig.  • The  new  approach 
has  already  allowed  for  an  extension  of  knowledge  shout  the 
temperature  distribution  and  heat  transfer  over  a very  wide 
range  of  the  Pr-numbor ; io'U  V * 

12k.  In  the  following  table  the  salient  features  of  the 
boundary  layer  theory  are  juxtaposed  with  regard  to  three 
milestones  of  its  evolution. 


1 

I 


■p*i 


1 


i 

I 


* 

\ 


i 

i 

i 

} 


I 


i 

i 


i 

i 

i 

l 

B' 

r 

u 


i 


SAVORD  Import  2959 


TABLE 


Tlissf  s £11;; tor.i * or  rrogtig  in  Bounfiftr^ 


lsL 


Main 

Period Initiators Gharao  tsrlsti.ee 


Haat  Transfer 
Problem  soluble  for 


Early : 

1%$ 


Prandtl, 

Taylor, 

eto. 


Dichotomy  of  th*  boundary 
layer  in  a purely  laminar 
region  near  the  wall  and  a 
full*  nii*hn}ant,  Oti  fcer  region. 
Prerequisites  of  this  theory 
are:  (1)  Similarity  of 

boundary  conditions  in  the 
equations  of  momentum  and  heat. 
(S)  Ko  pressure  gradient  (flow 
■long  flet  plats).  - Otherwise 
a fictitious  hsat  souroe 
arrangement  should  bs  super- 
posed.  - 


Pluid  with  Pr'  ■ 1» 
However,  origlmal- 
ly,  both  *****  and 
Pr  were  taken  »s 
unity.  With  the 
eotual  vaiui t 

> ! 

the  original 
results  held  for 
the  single  value i 

pP  « J-  <1 
~ oj/f  f*  '’***/V 


Kiddle « 

1938/ 

1949 


Von  Kerman,  Triohotomy  of  the  boundary 
Relohardt  layer  in  one  purely  laminar 

eto.  wall  region,  one  intermediate 

buffer  layer  ( simple  or  oom- 
posit*}  and  one  outer,  purely 
turbulent  region. 


Fluids  with;  Pr1 
► 1 and  Vr\*^  “ 1 
The  theory  ho Id a 
on  a rather  limited 
range  of  walnea  for 
the  Pr -number . 


I 

1 

i 


i 

i 


5 

1 


4 


l 


Aotuali 

alnoa 

1949 


Relohardt. 

Blaer, 

eto. 


Continual  decay  of  th*  sffaota  Fluids  with  rr*>»l 

of  turbulenoe  from  the  free  and  such  with 

stream  to  their  suppression  at  Pr'«<l,  and 

the  well,  A characteristic  variable  ■ on  th* 

width  ( *!,  ) , ia  re  fax  red  to  ths  ride  range  of; 

transport  of  momentum.  This  a 

width  is  defined  ei  that  dia-  **.  o®  ^ is 

tance  from  the  wail  within  - - 

which  Newtonian  shear  streea 

prevalla  ovar  tha  apparent  om 

from  turbulence, 

•Similarly,  for  the  thermal 
boundary  layer,  another  char- 
acteristic width  ( *}t  ) marks 
the  region  of  prevalence  for  the 
heat  transport  by  molecular  conduction 
over  tho  counterpart  due  to  turbulent 
oonveotion.  Except  for  Pr'  = 1,  the?* 
i*l  “*1  ki  tj. 

•i  * •» 

*2 
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Differential  Equation  for  tho  Temperature  Distribution 

125.  The  differential  equation  for  the  thermal  profile  near 

the  solid  boundary  may  assume  either  form  (A)  of  the 
temperature  gradient:  ±Hz!  **  y (yj 

iMf 

or  form  (B)  as:  ctfe*.)  ^ <+>/«.) 

'dLu. 

126.  In  the  form  (A)  the  distance  from  the  vail  is  chosen  as 

argument  of  the  sought  function  with  either  both  variables 
being  dimensional  in  ? (y)  or  both  being  reduced  to  dimension- 
less varieties  like  { \./n  ) or  xt  ( yaV*/  ).  The 

form  (A)  seems  preferable  for  fluids  with  ?r*  « a* 

127.  The  form  (B)  of  the  temperature  function  has  the  velocity 

aa  argument.  Again,  either  both  variables  are  taken  in  the 
dimensional  form  T W)  or  dimensionless  varieties  arc  used  ^or 
both  in  ^ { 7 ) op  ^ ^ ).  In  most  cases  of  practical 

interest:  v/y  m Mfa*  \o  to  if  (Reference  (m^) ) • 

12Q.  The  form  (B)  seems  preferable,  for  Pr»  » 1,  because  of 
the  rapid  decay  of  the  integrand  function  with  the  dimension- 
lees  distance  ( '*//%.  ) from  ths  wall. 


*.  Form  A 

129a  With  reference  to  formulas  (19),  she  total  flux  density 
for  the  heat  may  be  expressed  as  follows* 

o ^0  + *)  ^ 

which,  with:  a ^ A r*  W fi  , k,  *1  yields: 

n ( \ 4*  r ?k 1 • «*T 

I30.  By  accounting  for  tho  boundary  value  at  the  wall; 

a 1 t AT  \ and  markinf?  a symbolic  Integration - 

'•  ’ M *fl.  " ~ 

one  obtains : 


T-T 


r 7- 

_ [HI  i. 

\ A?  /.  1 


(27) 


1 • V»  i vil  d liuc  hoi  '."'.'o##  ViriBOiea  1 ^ »|T“V  ’:5?  ano  y/ r . the 

relation  Is  transformed  Into:  *‘w 


33 
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r W* 


ft-  / rt>  \ 

^ - I 2i) 

» u 

M 

fe. 

k 

0 

0 £ >/  n 4 l — 

^ r>  c 

~r  w — 

-j-  i » 

b. 

Form  B 

(28) 


132.  With  reference  to  formulas  (19)  the  expressions  for  0, 
and  may  ba  combined  into  j ' m 


AT  m £_ 


On  the  otner  hand: 


fc  T- 


(29) 


Q-l 
» ’£ 


And  -cl~ing  for 


* prz- 

fc, x 


133.  By  substituting  this  expression  for  , the  equation 
(29)  ill  transformed  into  t 

y _ ~~'r (30) 


dT  _ _K 

«Ut 


L — _ T?  1 


» ♦ Jl-  (.  ?-  * 0 

l-hf\ 


13l4.,  By  taking  in  account  the  boundary  conditions  at  the  wall 
( 't,  *.0  ^ ) * 

AX  JU  9, 

**  <«.  r. 

The  symbolic  integration  of  equation  (30)  yields: 


t-to  -fg  [ 

V v*  U 1 1 t l w - t 1 |""f  1 > ♦ /»  w- 

J.  M-A  * # 


(3D 


for 


r 


K L 


(yj)1 


JL 

V * 

■ • 
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135.  With  the  influential  properties  of  the  fluid  assumed 
constant-  1-e.  with  1 » I , and  a shift  to  dimensionless 
variables  ( v slT-T^iV »n<3  <f  * <u./V  ),  the  differential 
equation  (31)  is  reduced  to: 


(32) 


Integration  of  the  Temperature  Equation 

I36.  The  integration  of  the  temperature  profile  is  discussed 
for  the  form  (B)  of  the  differential  equation  in  its  dimension- 
less form  r if  ).  The  same  method,  however,  is  easily 
applicable  to  ^orm  (A)  as  well. 

137«  In  a slightly  modified  way,  equation  (3?)  may  read* 


v*  - (£U  lU  V 

I I*  l Jm  J*o  f J 


and  by  means  cf  the  identity: 

> v p w 

equation  (32a)  is  transformed  into; 


_ 1 


,NP*' 


(32a) 


(33) 


v * 


138.  The  important  function 
a Binomial : 


may  be  expressed  by 
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with  s as  a dimensionless  quantity  which  nearly  vanishes 
at  tbjs  wall  ( J\m  °),  while  assuming  a value  ©<A 
afar  from  the  wkIjl  . 


I39.  The  value  ^ = *.  * 1 may  be  used  as  a first  approxi- 
mation to  yield  ( )*  which,  close  to  the  wall,  la 

identical  to  the  function  T/- % « • The  latter  is  a linear 

function  of  the  distance  from  the  wall  for  the  prismatic  ducts 
that  are  considered! 

y,  - 1 - ik  - Vn.  (35) 

' 9 


1 i.n 


Tf 


Wf-  4. 


Wi  rs* 


ij  wiiv  iuuy  vavi!  ' # M.jy  D9  round  r^oss 
the  continuity  of  heat  flow  by  balancing  out  the  amounts  of  heat 
involved  in  flow-directions  parallel  to  the  mean  motion  and 
across.  This  reflects  a similarity  between  the  temperature 
profiles  at  different  stations  along  a cylindrical  duct.  ?or 
instance,  this  means  for  a circular  tube,  that,  over  the  radius, 
the  flux  density  of  heat  is  accounting  for  a proportionality 
between  amounts  of  heat  that  are  transferred  radially  and 
axially. 


ll±l.  Thus,  the  dimensionless  heat  flux  density  ( ) may- 

be expressed  by  the  following  integrals  (Reference  (m^) ) 1 


Tube  x 


Channel : 
(rectangular) 


/..I  - 

* !?/  1 


g'y  C » - ’/*  ) *K'V 


'll  \ 

/*•  I 


^ o 


< 


0-  M'*  I 

I 


(36.) 


I - 


a K/% 

J rf^ 


\ U¥M 

I I \ / 


1 


-J 


A( 


v/  V 
• #e  \ 
v } 


(36b) 
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HI 

* 


t 
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„ t ^ 


1 42.  In  replacement  of  * by 
tomperaturo  profile  read's  T/<« 


, the  equation  for  the 


Jr 


l 


dt) 


ft 

i f/i 


U<f/«  j 


f/t(>  A-  /*>  ft!  ) 


<kjf 


(37) 


Provided  the  velocity  profile  ( P ) and  shear  distribution 
{ 'r!ta  ) are  known  functions,  the  temperature  profile  ( ~7"  ) 

may  be  determined  from  the  equations  (36)  *nd  (37)  by  a method 
of  successive  approximations: 

II4.3.  As  first  approximation  ( ) for  the  thermal  field  in 

the  boundary  layer,  a functional  fora  may  be  chosen  to  correspond 
to  \ 


A » Vi 

"o 


IQ/  \ ft-  / 

\ } * -/T- 

'#  / I * 

* 


Thus  j 


<v  (H\  ( j_dL _<u 
\*rl.  j « >-[■*  n 


t 


" ['  " being  a known  function  of  1/\,  and  , while,  for  the 

unknown  function  Pr’  * mPr,  a convenient  average  value  over 
the  boundary  layer  may  be  estimated. 

lllij..  If  la  entered  in  one  of  the  equations  {36a)  or  (38b) 
together  with  the  known  velocity  distribution  C */*.  ) ■ a 

second  approximation  ( )*.  can  be  determined  regarding 

the  heat  flux  density,  a ay  ror  the  channel : 

•i  » 

7/X 


r 


la/  \ 

i % i. 


! Cv~  A { ?/*  \ 

J V,  f 


['  7,f  <u yM 


A second  approximation  ( ) for  the  thermal  field  will  be : 

Y . 

94*  \ ( : ±-  /*  1 

J - (33*) 


T 

V. 

A. 


\ 

/e  J ^/r.  C 1 


17 

«/  I 


I I 


xr  A 'trf.  TJTS  rs 


V?,:||+^  . 

2 1 <W 

A p < I 

\ + fin. 


(iiOb) 


Integral  In 


II4.9.  In  the  first  part  of  this  review,  the  dimensionless  , / . 

velocity  profile  has  been  expressed  as  y ) , with  V * "Vac.  • 

150 s Sy  transforming  the  integral  Ii  from  u to,  f , and 
using  the  identity: 


Q8L 

<h 


TL 

n9 

» t fi 


we  obtain: 


\ ■*/«  V V 


dU>  =,  — 


•T- 

-U.  J Hi 
CAT 


I . ,\f.  . m 0 t\ 

*j  0kv+r"v*‘rr,*-f 


4 / 

151«  As  the  function  ft  * ^13  sharply  growing  with  the  dimen- 
sionless distance  in  tho  Reynolds  form:  ->7  •»  1 •**/*/  , the 

integrand  is  declining  with  growing  . Hence  I,  $ 0 only 

very  near  to  the  wall  where  r « T,  . This  is  expressed  in 

"I  4 4.  J _ i-V,  ~ rrM . 1 T f U ^ fl  •» 

Ci.USSOi.Jf  ij  Aiwiu  A4«^  V.UO  Xii  uo  9-A  A 1 »»  •*.  WA4  —•jj  * 


! 


i n 


i * 
\ » 
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»>*  Integral  1^ 


I 


v ?J  r t + a ) 


duf 


lf>2.  Very  near  to  the  wall: 


0 


1 • e • 


I,  C 

4® 


Afar  from  the  wall;  / - { j hence  the  unity  bs;comea 


negligible  as  compared  to  /*  and  to  ^ ft-  (for  Pr'  not  toe 
small),  and  for  high  values  of  y (or  one  of  its  dimensionless 
varieties  we  obtain* 


**(•*/>> 


/V 


l 4-  fiK 


m 


( 1.  a \ 


S - 


« I 


« t 


’9 


relevant  founds.: 


j ccnd  itiona  it  may 


153*  By  considering  the 

be  shown  that  the  integral:  v - j > -f  . over  the  whole 

boundary  layer  would  be  negligibly  small  for  a flow  past  a 
flat  plate.  The  same  integral  would  b*  le--  negligible  for  a 
two  dimensional  flow  in  an  open  duct,  and  still  larger  for  the 
flow  in  a tube.  This  means: 


t 


'pin  UkMUi( 


r..  /.+ 


Integral  Form  of  the  Temperature  Profile 


154*  With  the  expressions  ^ and  , the  integral  form  for 

the  temperature  profile  ( ii.0 ) becomes:  ’ 


S-  , 


4 Mostly  the  symbolism"  of  Reference  (m)  is  adopted.  The 
exceptions  include  rases  where  a conflict  had  to  be  prevented 
with  well  accepted  American  practice,  e.g.  by  using  r,kn 


Lnstesd  of  for  the  hest  ecn 


i du  c uivx  t»y . 


I r» 


r? 


HAVOHD  Report  2958 


e»  1 a.  I dl5  \ r «« 

»T.  ^ « i ~£~  | L.  r 

\ lo 


_ (fJ-i)  a.  •*■ 


n 

V.I 

f J 


(W*) 


155*  In  da no ting  simply  by  "an  and  w (without  subscript) 
tha  raspective  values  of  the  two  integrals  at  a limit  distance 
from  the  wall  for  which  ( y *=  r)  t * y>  ».  I . vo  obtain  as 
gradient  of  the  temperature  profile  at  the  wall t 


[^L\ 

i, 


145) 


. + *iZ.  Cf-r'-O  » V ' 

1 * **  r 

v 


156.  Finally  the  Integral  form  of  the  temperature  profile 
reads : 


Cl- 

*0  » 


& I 

H-  (.fc'-i) 
<f  + \> 1 


(1*6) 


jlT  c?;  -0  *.  * 1 

1 u 


157*  Similarly,  with  \ii  : 

» • *^r 


#»  / 

4 + eP-'  -0  ~y  + /r 


'*  il'-AiL*.  * l 

■ ' - 'V~ 

Characteristic  plots  concerning  the  thermal  boundary  layer, 
as  obtained  by  integration,  are  shown  in  the  Figures  8 - 11 
(References  (m_ } and  (ny)).  The  following  indications  ;.rs  to 
supplement  their  captions: 


/ 1. 4_  \ 


1 
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Figure  8 These  profiles  are  a second  approximation,  ( £/* 

of  the  heat  flux  density.  * 

The  dash-dotted  parabola  represents  the  heat  flux 
density  for  the  laminar  flow,  a;  obtained  by  a 
similar  numerical  integration  in  two  steps. 

Figure  9 The  first  aDproximation  (dashed  line*)  re.fwT’s  to  a 
flat  plate,  the  second  one  (solid  lines)  to  a 
channel  duct. 

The  .small  rings  mark  the  points  whose  abscissae 
measure  the  width  M ''It  w of  the  prevalently  heat 
conducting  layer.  The  single  abscissae  " % M - 
independent  of  Pr*  - is  a limit  of  prevalence  for 
the  molecular  viscosity. 
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Fart  III:  The  Turbulent  Heat  Transfer  at 

Smooth  "Walla 


Introduction! 


158.  With  a geometry  of  flew  as  considered  here  (prismatic 
duct  or  plate),  the  analytic  expressions  for  the  function 
Ar!p.  - /*  (. 1 i over  the  several  regions  of  the  boundary- 
layer  - as  recently  discovered  (References  (m),  (q)  and  ( t ) ) - 
have  allowed  to  determine  more  rationally  the  profile  of  mean 
values  for  both  velocity  and  temperature.  The  distributions 
may  be  approximated  analytically  with  a local,  differential 
formulation  by  integrating  hereafter  either  directly  for  the 
velocity  (Part  I of  this  survey)  or-  stepwise  for  the  temperature 
profile  (Part  II). 

159*  A -chief  objective  at  present  13  to  discuss  the  heat 
transfer  problem  for  turbulent  flow  as  handled  by  means  of  an 
iterative  integration. 

160.  The  theoretical  results  from  a local  approach  are  compared 
with  those  of  global , semi- empirical  relations  for  the  heat 
transfer  coefficient  at  the  wall.  Ths  latter  appear  as  a 
product  of  powered  similarity  parameters  as  factors. 

161.  With  moderate  restrictions,  as  specified  in  References 
(m)  and  (t),  two  integrating  varieties  of  the  method  are 
evolved : 


a.  The  computing  or  iio  at  transfer  i3  based  on  detennlning 

first  the  temperature  distribution  over  the  boundary  layer 
(References  (m) , (q),  (s)  and  (fc))  by  successive  steps  of 
integration;  then,  the  slope  of  the  thermal  profile  at  the  wall 
is  computed  and  entered  into  any  one  of  several  accepted  for-n  s 
of  the  heat  transfer  coefficient,  e.g.  the  dimensional  one  (h) 
or  the  dimensionless  varieties  ( p/. ; Ku,  5t,  etc.). 

4.^ 

b.  The  heat  transfer  at  the  wall  may  bo  determined 
directly  from  a general  relation  for  the  temperature  gradient 
at  the  wall,  without  full  knowledge  of  an  integrated  tempera- 
ture profile. 

162.  The  first  method  is  based  on  evaluating  all  of  the 
integrals  involved  in  the  temperature  profile.  For  a two  step 
approximation,  their  number  is  four.  Two  of  them  ore  indefinite 
integrals,  and  two  are  definite. 


to 


*3 
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I6I4..  For  then,  the  limits  of  integration  are  at  both  and*  of 
the  boundary  film.  Integral*  with  a variable  upper  limit  would 
be  irrelevant  here,  although  they  are  indispensable  for  tracing 
the  complete  temperature  profile  from  point  to  point  along  the 
normal  to  the  wall. 

l65*  The si  method*  of  successive  integration  are  of  praotloal 


interest  ee  long  as  a so-called  universal  law  holds  for  the 
function  ft  ( b )*  This  would  mean  an  indspendenoy  of  /* 
from  the  individual  flow  conditions  in  regard  te  the  geometry 
cf  wall  configuration  and  properties  of  th®  fluid.  A 
restriction  in  this  sense  i*  associated  with  ths  Poclet  number 
not  being  allowed  to  underatep  a certain  lower  limit,  vis. 


166.  This  limitation  corresponds  to  flow  conditions  which 
coubine  in  such  a way  as  to  keep  the  characteristic  width 
* \ " of  the  prevalently  heat  conducting  layer  within  the 
range  of  validity  of  those  forms  (A  and  B)  of  the  funotion  ft 
( -n  ) tlat  hold  near  the  wall.  Beyond  this  range,  i.e.  for 

"*1^  > so  (References  (a)  and  (t)).  " “ would  depend 

not  only  upon  the  wall  distance  (say:  •»!  in  dimensionless  form), 
but  upon  the  R«-number  as  well,  through  the  quantity: 

**}  « (iV-m.  that  is  involved  in  the  C-form  of  the  function 

/i  (References  (m^-)  and  (t  )). 

167.  A check  on  whether  or  not  the  actual  Rs  is  high  enough 

1 for  a practical  validity  of  these  methods  of  integration  would 

Boom  especially  desirable  for  fluids  with  a very  small  Pr-number, 

! which  yield  high  values  cf  the  ratio  2500 /pr. 

li 

168.  On  the  other  hand,  for  a very  high  Pr-number,  the  integral* 
involved  in  the  heat  transfer  calculations  show  a trend  of  a 
considerable  simplification.  The  numerical  value*  of  t-h*  heat 
transfer  coefficient,  as  obtained  by  successive  integration*. 

j - were  found  In  good  agreement  i?ith  those  baaed  on  semi-empirioal 

| power  relations  as  reported  in  Reference  (m)  for-  tie  higk 

Pr-number  a *nd  in  Reference  (q  ) for  the  low  ?r-numbers.  Thus 
I * the  oompax-ed  results  cover  a wide  range  of  the  generalised 

li  Pr-number. 
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Analytic  Approach  by  Iterative  Integration 

I69.  The  profile  of  the  dimensionless  me  an  temperature  - 
say  ( >/*- ) of  a circular  tube  with  radius  r - may  be 

determined  for  nearly  constant  average  properties  of  a turbulent 
fluid,  by  successive  steps  of  integrating  the  following  two 
equations  s w 


Ilf  Cl  - Vt ) 


mwrw 

**1*  - VoaCA) 


r oyc-^)^a) 


— (U-7) 


[AL\ 

W. 


jl  /\  ) 

— — *y 

M i \ 

+ /i  r*: ; 


170.  two  unknown  functions!  and  T are  thus  inter- 

related to  each  other  and  to  three  nearly  known  expressions : 

/e  * '’f/n  , =-  **/lr  and  T/T.  • A dimensionless  either 

wall  distance  'll*.  « *1  * l^/>) or  velocity  ( y ) 

is  a common  argument  for  *11  of  the  five  functions  involved  in 
the  above  integral  equations.  As  to  Yt! 

it  is  a parameter  for  which  a mean  value  ov*r  the  boundary 
layer  is  used  according  to  the  2 s cmc  t»ry  of  tillS  flow  ( ^ ) 

and  the  nature  of  the  fluid  ( 

171.  As  first  step  or  approximation,,  the  expression  %*Yr. 

la  substituted  for  in  equation  ( Lpo ) which  then_  may 

yield  by  a single  integration,  the  first  approximation,  , 

for  the  temperature  profile. 

172.  The  second  step  consists  in  substituting  \7j  for  % in- 
to the  two  integrands  of  equation  (47)>  thus  obtaining 

after  having  evaluated  the  two  integrals.  The  second  approxi- 
mation of  the  heat  flux  density  is  used  in  turn  for  integrating 
equation  (I4.8)  once  more.  This  yields  the  second  approximation 
( ) or  tha  temperature  distribution. 


45 
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17 3*  iiy  substituting  iteratively  the  aucces8iv#  expressions 
for  { 5V9*  );and  ^2  from  one  into  the  other  integral 

equation,  the  step  of  approximation  for  the  temperature 
profile  may  be  reached.  It  would  imply  the  evaluation  of 
1 + 3 ( 'n.  - 1 ) integrals  involved  in  the  general  expression: 


‘ l*/Q> 

. . evo  (.  •*  /**•'> 

V 


(49! 


I7lf.  Except  in  casea  of  very  low  Pr-numbers  (molten  metals 
and  alleys,  etc),  two  steps  of  approximation.  i.e.  n ■=  2, 
should  suffice.  Then,  the  number  of  integrals  to  bs  evaluated 
would  be  1 + 3(2  — — — 1 ) s 1^. . 


175*  Ae  to  the  choice  of  the  dimensionless  argument,  either  a 
relative  velocity  ( «j > »**/V  ; * **■/*<  * ) or  a 

relative  distance  from  the  wall  ( 17  ■*.  m *V  j for  the 

circulni-  tube:  '•.in.  ) have  been  used  according  to  conditions 
of  flow. 


176.  With 


t 


! 

* + ( At/u  ) ^ icJl 

V M , 


two  out  of  four  integrals,  for  n *=  2,  have  a variable  upper 
limit  within  the  boundary  layer,  while  the  extent  of  integration, 
for  the  other  two,  reaches  the  free- stream  value  of  the 
relevant  argument. 


177.  The  two  indefinite  integrals  may  show  one  of  tha  following 

i*  o *rwTW  «s  • 


- f* 

O 


n 


(50) 


v5l> 


NAVGRD  Report  295° 

178#  Th«  corresponding  definite  integrals  with  subscrJ.pt 
denoting  an  extension  of  the  integration  over  the  entire 
boundary  layer  are j 


B 

Vtu 


A4r 


w«  * »**.  " • 


K 4^) 


fH(yp*Ml 

1 


(52) 


/ . . 

1 » 1 
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It 

1 w 

r«,  .tKt^  • 
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■T  jw  -r  dora) 


j (y/*i 

■'  I 


>d«i 

TT^ 


f C*  t 

a *.'«✓  / 


1'9 • The  expression  for  the  heat  flux  density  at  the  wall. 


(54) 


muy  be  transformed,  by  the  uae  of  dimensionless  variables: 

* 


•V  T-T_  «*/..  a*  - *1***  V » #AC  into 

m ~~=  , t V 'W  ' ''  " ” ’ • 


(55) 


100.  As  shown  in  the  second  part  of  this  review,  the  slope  of 
the  thermal  profile  s.t  the  wall  may  be  given  the  fora: 


71 
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l8l.  Similarly  it  was  shown  that,  in  a two  step  app.roximationp 
the  expression  for  the  profile  of  the  moan  dimensionless 
temperature  ( J ) is  finally: 


(%). 


p * (Rl-.)I  ♦ U „ ■ (57) 

I 


182.  By  replacing  in  relation  (55)#  the  slope 
by  the  expression  (58),  one  obtains: 


feG  -r  ~ -P. 

v m 


<i  * sz — — ."r: — 

C + {Psl  — i)l  ♦ 

<t  ** 


or  \7 0 vo  \v*_ 

V > • 

' ♦ ^ ')L  * te 


(59) 


183.  The  accepted  forms  for  the  heat  transfer  coefficient . 
at  the  wall  are  either  dimensional  ( K.  without  • superscript) 

or  dimensionl-sflR  ax#  1*  *■  *K.a  /#  * 

* ^ * - ({• 

Cf  rf.  *V  . ,4.'  _ i!» 

St  xTK  ■«  > * ~ xX  nr 

I04.  The  subscript  t refers  to  the  temperature  difference 
\ iJ'i  ) used  as  a divisor  An  deriving  the  dimensional  heat 
transfer  coefficient  fro®  the  heat  flux  density  at  the  wall 
>.  0 )»  The  second  subscript  ( j ) Indicates  the  reference 

velocity.  Below  are  listed  seme  of  the  current  varieties  in 
forming  the  dimensional  ( k; ) and  the  dimensionless  form  ( ) 

cf  the  heat  transfer  coefficient  at  the  wall:  ^ 


Dimensional  form: 


u 

1 *T- 

1 


h.  * t h-x. 

Hk  % I T,-T. 

h , % iX 

VII 


(59) 


* 


Dimensionless  form: 
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K % 


*Lir 


**!,V 


a #>  1 T if1 * * 

' V " 


(60) 


185.  Two  varieties  are  used  for  the  second  subscript  n J-  * 
by  taking  its  reference  velocity  either  the  free-atrsam  value 
( V ) or  the  shear  stress  velocity  ( ^ viz: 


% 


‘f*4  !Ah 


* 1 


’ a c 0 V 
<d\J  »« p, 

anu  similarly: 


» I \ ^ *>*»*/• 

•g  '* 


(61) 


V, 


UiVv 


ijfl 

IT" 


= — !.'(£)  - (62) 

•>f*  - «V  ’ *r  =■* 


with 


1C 


p*  4^1 


f*  r*  fit1-.  J_  ^ 

if -1  _ I i<  - 1 - 1 fit  dl*At 

?*'-«  ' J *<,7'  Jo**  1 

, M ( O ^ w 


V 


I 

186.  Except  V>;/  » ^wo  ®ore  dimensionless  varieties . are  in 
use  for  the  heat  transfer  coefficient,  vie.  the  Nusselt 

number  (Nu)  and  the  Stun ton  number  (St).  These  may  be  expressed 

as  follows: 


Nw. 


v * 


i. 


Ha 

' Hu 

^ - 


^ ad 


( &*) 


^ * - fV  ' P.crU~tT"  ' ¥„*«.  *»1,‘ 


Pi 


U-9 


H 


<» 
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I87®  A general  formula  for  the  Nu-number  may  be  derived  by 
combining  the  following  two  equations : 

it*  xfLt. 


u. 

*3r  , 

*•  f iii 

f iL  \ | is. 

\AfW/„  J.  K 

4yA) 

1 + Ar*. 

i 

^ ^ S*  1 

▲ rcm  wiis  bcundsrj  condi cion  ^ ^ for 

temperature  gradient  at  the  wall  is: 

i U \ 

1 

V <wi. 

V ^ hit 
J.  T *f 

- A 07*) 

R.' 

138.  On  the  other  hand,  with 

* s e I \ 

tsV  « TA  Qm  - K\T7 

m.  ~ ~ *e  \~7  l 

It  follows  that: 


Jr. 


\&) 


r *•• 

*!*•  <k 

1 fe 

«««' 

£614.) 


189®  The  main  integral  in  the  denominator  of  formula  (62)  may 


w-.  - — _ ^ _ 4.  . j uM  . 

uo  a|.’yi'U Axiua  u^u  ; 


! T \ 
[ \ ! 


r 4.- 


r . 


\ C m 


Ay/A ) 


190.  Thus  for  a logarithmic  profile  of  velocity,  Implying 
<kg  o » integral  becomes: 

r X (YL_\ 

~vf 


I 
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191*  It  is  interesting  to  note  that  the  integrand  of  ( A,,} 

* approx. 

presents  a more  rapid  decav  with  the  dimanaionlea*  wall 
distance  than  would  occur  with  the  Integrand  of  the  denominator 
f ‘ JJSx  formula  ( 64) . 

\ ' 

192.  Consequently  with  the  aid  of  the  relations  (6l)  or  iu2) 
it  may  surfice  to  apply  the  "A“  and  nB"  form  for  the  function. 

l vj  ),  provided  the  Pr-nuraber  is  not  too  lew. 

193.  The  general  expressions  (formulas  (6l).  (62)  and  (63))for 

the  dimensionless  heat  transfer  coefficient  ( AW  ) »ay  be 
specialised  for  the  conditions  of  the  early  theorv.  It  implies 
Prandtl's  dichotomy  of  the  boundary  layer  in  a purely  laminar 
wall  region  ( /*v  - 0 ) and  a purely  turbulent  outer  region  , T 
(a  •-=  ).  Thus,  with  the  particular  values  *!.•>*  *,**'tS« 

*'**4/'  the  relations  (6i)  and  (52)  of  the  advanced  tfTeciry  would 
become  t 


U ) 


( 61  * ) 


i 9 


■ I#  • • ■ n 

\ + \jt  ~ • / T* 


U \ 

V Vi 


t**jL  •H \ -V  ~ *f( 


(52*) 


19m.*  In  these  formulas,  which  denotes  the  relative 

velocity  &v  the  border  between  the  two  regions,  is  a purely 
hydrodynamical  quantity.  By  comparing  the  particular  formulas 
(primed)  to  the 'general  relations  (6l)  and  (62)  it  follows 
that,  in  the  actual  third  phase  of  the  boundary  layer  theory 
for  heat  transfer,  the  counterpart  of  v#  ia  the  Integral! 


* f &'-1  A* 

f i K-'  "ft' 

V 


f'  P-i 

J Rr'-  % 


195*  Evan  by  simplifying  with  a constant  ( 1 ) t inataad 

of  including  the  variable  ( ft/  — . i ) 5.n  the  integrand,  still 
the  simplified  integral*  M f*  will  remain 

'f  ) HTpT'  * J !*aR’ 
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a mixed  thermo-aerodynamic  quantity.  Tills  is  due  to  its 
dependency  in  the  denominator  upon  the  generalised. Pr and tl~ 
number  ( 'VJ  ) multiplied  with  the  function  * 

Physically,  the  integral  f* dy/i«7*fr  expresses  the  relative 

i i 

velocity  for  -a  » « i;e.  et  the  relative  distance 

from" the  wall  for  which  ^ ^ . There,  heat  motion  by 

molecular  agit&t  ion  and  by  zncl&r  convection  due  to  turbulence 
contribute  equally  to  the  heat  transfer. 

196.  Prom  the  above  analysis  two  method®  of  successive 
integration  are  evolved * 

a.  The  heat  transfer  calculations  consist  in  determining 
first-  by  iterative  steps  of  integration  the  distribution  of 
both  the  temperature  (function  T , formula  (4-9)  OI*  (5? I)  and 
its  gradient  over  the  boundary  layer.  The  wall  value  of  the 

tup#  gradient  is  then  eomp'iatad  and  entered  In  the 
relation  for  th®  heat  transfer  coefficient  (formulas  (56)  to 
(63),  References  (o),  (m) , (r),  (■))• 

b.  The  heat  transfer  at  the  wall  may  be  determined 
directly  from  the  general  relation  for-  the  temperature  gradient 
at  the  wall  {formula  (5&))  without  soaking  previously  » full 
knowledge  of  an  integrated  thermal  profile. 


197.  The  first  method  l~  based  on  evaluating  all  of  the 
integrals  involved  in  th#  temperature  profile.  For  a two  stsp 
approximation  their  number  is  four.  Two  of  them  (Ij,  I2)  a**s 
indefinite  integrals,  and  two  are  definite  (1^  , I2  / • 


I98  The  second  method  implies  the  evaluation  of  the  two 
definite  integrals  only.  These  are  (formulas  (62))  I]_|  and 
lo*  , end  may  be  simplified  by  taking  the  factor  iv-i  outside 
the  integral  as  — I . The  limits  of  integration  are  at 

both  ends  of  the  thermal  boundary  layer*  Integral#  with  a 
variable  upper  Xisit  are  irrelevant  here,  ilthoush  they  should 
be  indispensable  for  tracing  the  complete  profile  of  the  mean 
temperature  from  point  to  point  along  the  normal ^ to  the  wall. 

199.  As  for  the  integral;  ^-l£  * 1 x * f 

? 


T/r9 


is  stronglv  deoendsnt  uoon  the  geometrical 


configuration  of  the  well. 


200.  For  the  flat  plate:  f «l  and  8-*o#  for  the  open  duct: 
fc  though  small,  is  not  negligible,  for  the  tube:  J y 1 einoe, 
at  a closing  wall,  the  gradients  of  sheer  stress  and  heat  flow 
density  are  of  opposite  sign,  and  the  relevant  profiles  of 
f ( ^ ) and  <j,(  7 ) may  differ  considerably. 
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Thus, 


vhe  most  important  of* 

sr 


\ 


it 


~ T" 
1 


j V -7 


T“  Vi  A •?  — _ *1  - * . — _•»  # _ _ 

ni.3  iiivuiVtiu  13  : 


201.  Afl  the  function 
poorly  known  ratio  ) f 

the  time  being  by  the  avers 
confined  to  the  integral: 


J?*aA  Wk 


it  might  as  i*6 11 


depends  on  the 
be  replaced  for 
Sad  tiio  discussion 


A 


.M/ 


t / 4-f  \ 


l ^2) 


202.  The  form  of  the  integrand  is  dependent  upon  the  width 
( ”7*_  ) of  the  prevalently  heat  conducting  iaver  near  the  wall. 

" Tt  " ia  that  characteristic  value  of  the  variable  for  which 
the  function  /*  ( ^ ) equals  the  reciprocal  of  the  preneralized 
Pr and tl -number : 

*>)«.)  '**-  */  Pe!  (66) 


203,  As  rrss  mentioned  in  Part  II  of  this  review  (Reference 
(ny*)),  the  function  /S  is  independent  of  the  Re-number  ns 
long  as  ^ < 30  , and  in  the  range  j fe  < >,  <.  3.  ? , the  form  of 

(*  is  : 


r> 

1 


(6?) 


whore  9 , while 

viscous  boundary  layer  i.e. 

20I4..  For  a region  very  close 


is  a measure  of  the 


to  the  wall : 


width  for  the 


lim  tanh  6 
\t  — > v-ery  small 
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while  afar  from  the  wall; 

lisi  tanh  S'  1 
G ->  very  large 

205-  For  small  Pr-mmbers , " Ij. n expands  into  the  region  where 
the  exchange  of  momentum  la  almost  of  a purely  turbulent 
nature,  and,  with  ■»>  and  6* «■  V*.  tend  ing  to  grow  very  large, 
so  does  the  Re-number.  Then: 


(■*)  - ^ "i 


very  large 


206.  Solving  now  the  equation:  V((^  — ^ r '/ftj  yields; 


/ \ 

V-  \ i -h 


VCp» 


and  the  relative  velocity  function  afar  from  the  wall  becomes: 


1 La.  *1 

- **  I 


‘ fc 


+ v. 


207.  With  the  definition: 


T * ( ^ )a  ' f K X! 

'iv  1 /*»i 


I I A.  11.  II  \ I J 

'■K  1 />»| 

the  validity  of  relation  (7C)  is  limited  to  regions  with  a high 
enou  gh  value  of  tne  Re-number. 

208.  For  formulas  (oil)# 


l,  ‘j 

At  ~ V* 

u ^ 


ft- 


and  with 


...  iv  >1 

" “ "*  'h»  1/ 


it  follows  that; 


' ZY  >» 

"N  *n  ^ 

209 « For  a small  Pr,  the  term  VkR'  grows  so  large  that, 
comparatively,  the  constant  may  be  neglected.  Sene© ; 


V 

vi.) 


t / .s' 
IK  n 


Iwiti  ft 
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It  follows 
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210.  Close  to  che  wall*  ‘vw  ~ ' i.e.  ?f  ~ If 
that  for  3iuall  Pr  and  large  Re: 


v — csf 


^ 1 <u 

i nn 


r>  _ 

itn 

ln. 


(75) 


K 


211,  It  was  found  under  the  conditions  envisaged  in  Reference 
(m  ),  that  the  validity  limit  of  the  B-form  (formula  (67)) 
corresponds  to  7/*.  •x  ° r . This  figure  is  based  on  the 
best-  emp^ical  data  on  turbulent  flow  in  trie  tube.  Hence, 
with  ■—  o.  S’  and  K » u.S  , the  method  of  direct 

integration  for  the  heat  transfer  coefficient  at  the  wall 
would  be  limited  to : 


•1 


(■ 
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Vtc^V}  ) 4 
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212.  It  follows  that,  with 
Vi*_**U-  > and  with 


«■ 


\ - V 

V ~ if 


CA 


r<fc 


(Reference  (nij)  ) 


TU  ^ 00 


{76) 


SemVempirical  Approach  by  Power  Relations 

213.  Between  differing  conditions  of  a flow  with  heat  transfer, 
the  general  criterion  of  similarity,  as  based  on  dimensional 
analysis,  appears  as  a function  " by  which  the  heat  transfer 
coefficient  at  the  wall  (in  any  form,  say:  Nu ) is  Interrelated 

with  six  parameters,  vis: 


N,u_ 


{ ( 1. , 1 , s, ,K ,K,K 


1 77 ) 


214.  The  arguments  are : 

P . 1?.  _ 


: Reynolds  number 


Generalized  Pranatl-number 


Viscosity  ratio  with  subscript  fl 
referring  to  the  msitn  temperature 
oj  uiic  x J.uid  Hiid  w vo  tne  mix  • 


t QTMlOt  r»  f"  4 /sf*  f Ha  /v4 

- • -v-  W .1  w v-v/  v/*  Wtfv  '•»  ».*  • - V 

closed) 


(open  o v 


Recovery  ratio  between  the  equilib- 
rium temperature  iT.)  at  the  wall 
and  its  stagnation  temperature  (Ts^). 


215.  In  time  prst  it  was  observed  that,  oeteris  paribus,  th® 
heat  transfer  with  streaming  gases,  if  compared  to  that  with 
viscous  liquids,  appeared  inversely  effected  by  a reversal  of 
the  heat  flow  direction.  Those  experimental  facts  were  rightly 
correlated  to  the  opposition  of  sign  as  for  temperature 
coefficients  of  viscosity  between  gases  and  liqiida.  Independ- 
ency of  the  r:>te  of  heat  transfer  from  iv.a  direction  was 
secured  by  intu  during  a new  argument  (i\)  and  txus  correlating, 
ceteris  paribus.,  vi  “o  a power  of  P with  *n  empirically  found 
exponent . 

216.  The  parameter  * ! 4^1  accounts  for  a variable 

viscosity  over  the  temperature  profile.  By  neglecting  P„ . the 
early  investigations  became  affected  with  a dependency  ofths 
best  transfer  coefficient  at  the  wall  upon  the  flow  direction 
oi  ins  ixOau  1 1 w s upon  whether  hes<*  was  moving  toward  the  wall 
or  s. way  from  it. 


217 • The  parameter  = x«/Tstagn  be  implicitly  introduced 
by  taking  the  difference  (T9  - Tq)  between  the  equilibrium 
temperature  ( T0  ) at  the  wail  and  its  actual  temperature  (T0) 
as  reference  in  forming  the  heat  transfer  coefficient.  Them 


t ('Re  , Pfcl,  ^Vk.  ,*/l-  . H ) 


{?S) 


f 


NAVORD  Report  2958 

218.  The  Mach  number  extends  the  similarity  requirements 
over  the  compared  density  fields  provided  the  ratio  ef*/c»r 
linked  to  the  Pr -number  by  ■ well  defined  relation  as  it 
occurs  for  the  gasss  at  least.  For  lower  speed  ranges  the 
dependency  on  the  Mach  number  he  neglected. 


219.  The  semi -empirical  prooedure  consists  in  giving  to  the 
unknown  similarity  function  “ 1 “ the  form  of  a product  of 
powered  individual  parameters,  and  trying  to  determine 
empirically  tile  unknown  exponents.  As  a rule,  power  relations 
with  a rather  restricted  range  of  validity  are  linked  to  e 
number  of  ( 1 + K ) empirical  constant*  with  15  N n denoting 
the  number  of  independent  similarity  parameters  ( P.  ) 
involved.  Then  Nu  appears  in  the  form:  *■ 


tu  . c.  ?;•  c p;*. 


(79) 


220.  In  order  to  reduce  the  restrictions  about  the  range  of 

validity  of  such  power  relations  it  had  been  attempted  to  cj- 
generalize  by  using  binomial  expressions  of  the  form  ^ ^ 

instead  of  £.**■  , Then  the  number  of  empirical  constants  In 

the  generalized  power  relation  exceeds  that  of  the  simple  power 
formula  by  the  number  of  factors  in  binomial  form. 

221.  Of  course  increasing  the  number  of  empirical  constants 
does  not  eas«  in  general  the  work  of  determining  them  from  a 
group  of  experimental  series.  Yet,  for-  independency  of  the 
heat  transfer  from  the  heat  flow  direction,  ceteris  paribus, 
still  the  following  generalized  power  relation  (Reference  (o)) 
would  appear  as  one  of  the  best  available: 


rJii  _ «'  h£ 


■*y»  _ y>  —» o It  r ni*  1 

r f -iLr\?  r 0 ♦ \ 1 

>•  1 • / *.  > 


222,  The  parameters  ?i,  rg,  P^.  are 
between  the  wall  temperature  ( Tw  ) 


th*  fluid  { m ^ 

..  _AJ_  . 


referred  to  the  mean  value 
and  the  mean  temperature  of 


22 ji.  The  coefficient  Q.11&  has  boon  50  chosen  that,  for  an 
aspect  ratio  of  D/L  = 1/150,  equation  (SO)  is  reduced  to  the 
form  of  a plane  mean  flow.  Then: 
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ZiU±»  A co mao n weakness  of  all  the  3 arh --empirical  heat  transfer 
relations  5.s  the  limited  and  approximate  validity  of  such 
formulas  within  rather  restricted  ranges  of  the  variables 
involved 

225«  Refsrenoe  (r)  offers  a juxtaposition  with  a discussion 
about  ranges  of  validity  for  a great  number  of  more  or  loss 
empirical  formula#’  of  hast  transfer  cos-sutiti ons  in  the  fcmn 
of  a Nu  number  or  a Stanion-number : St  » fin.  i d! 

2^tb«  The  asune  reference  presents  «il dG ■ for  turbulent  hast 
transfer  in  the  tube,  a useful  three-seals  nomogram  by  which 
the  S t— number  * '&fr  (82)  i«  correlated  w the  Re- 

number and  the  Pr—  number.  •«#  w. 


d 

h 

T1 

*2 

T 

«n 


diameter 

length 

inlet  temperature 
outlet  temperature 
mean  temperature 
wall  temperature 


a.  of  fluid 


227.  In  this  nomogram  the  three  variables  cover  respectively 
the  following  ranges: 
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Comparison  be two 

vu 

Iwal  Wr  Ti  x p 

Approach  and 

228.  To  this  end,  th*  defirite  integral:  a = ^ &V 


has  been  evaluated  (hefersr.ee  (my))  over  the  range  ,1°  a ~ 
with  seven  different  value  3 of  the  parameter  Pr  Uo>u',\iT*  t,  io_  i©*',  ) 

For  each  one,  the  function  LV)  was  plotted.  The 

integrals  were  determined  by  measuring  out  the  relevant  areas. 


a.  Barge  of  Median  and  Higher  Pr-Number 

2^9 • In  regions  in  wt* ich  ths  integral  nan  is  independent  of 
the  Re-number,  the  relevant  theoretical  resul  vu  were  approxx* 
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b:  K^aRH 


NAVORD  Keror*-  ?qch 

mated  (Reference  (e  ))  by  n po®»r  relation;  athsor  °“ 
Pr-nuiaber  alone,  according  to  the  table  below. 


r 

To  tai  range 

of  Pr 

r 

rr-regions 

Range  of 

m 

Theoretical 
Integral j 

•theor. 

Semi  = 
empirical 
relation  of 
tan  Bo  5 eh 

median  Values 

/ ...  * m 

a y 01 

/.»x  - -0.30 

V jtu } r? 

-i.Jf 

Q,  e »)  ft 

-J.  ,-f. 

1 ^ Pp  i 200 

high  value* 

>»  <6 

(9.12)Pr“0*20 

230i  Tims,  the  theoretical  power  relations  in  ths  regions  of 
X'alldity  of  the  B-fonn  or  the  A-form  of  function  /*  respectively 
(References  (a)  and  (t  })  appear  in  fair  agreement  with  ths 
corresponding  semi-empirical  formula  with  two  constants  as 
determined  by  Ten  Bosch  (Reference  (m^-)). 


b.  Range  of  Lower  Pr-  Numbers 

231.  As  to  measurements  of  turbulent  heat  transfer  for  very 
low  Pr-numbera,  the  earlier  experiments  havo  been  gravely 
handicapped  because  of  the  unsatisfactory  wotting  of  the  solid 
surface  by  auch  fluids  as  mercury' flowing  over  current  wall- 
materials.  Those  conditions  yielded  much  t»o  low  heat  transfer 
coefficients. 


232«  On  the  other  hand,  more  recent  measurements  (Reference 
(s))  achieved  notable  progress  by  a more  suitable  choice  of 
substances  for  the  fluid  and  the  wall  respectively.  Thus,  the 
low  melting  alloy  of  sodium  and  potassium  in  equal  parts 
ensured  a full  wetting  over  a rail  of  pure  nickel. 

2j3 » In  Reference  (g).  the  comparison  between  experiment  si 

result*  * >d  the  thaorstl  (lo  [ f* nTwnnl  a 


o.% 

tU  - 7 f \0.*iS)  fi 


(82) 
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revealed  a satisfactory  agreement  for  ^ s Below  this 

upper  limit,  however s the  measured  values  appeared  as  much  as 
30£  under  the  computed  ones. 

234.  Nevertheless^  in  Reference  (q.*)*  the  comparative  discussion 
of  the  results  01  Reference  (s)  tss  rsna jred  in  the  light  of  a 
generalized  Pr andfcl- mimb«r - Thus  with  g±  * Re  m Pr, 
the  formula  for  the  turbulent  heat  transfer  he came t 

o.?  . 

Mvc  m V (v«.  OtI  >,  fV  *°3) 

iwAl. 

235*  It  was  shorn  there  that  this  formula  is  well  satisfied 
by  both  the  theoretical  results  from  iterative  integration  and 
the  experimental  ones  within  the  rang#  of  tha  measurements 
(Reference  (s)). 


WATrnon  f Ortn/i 

- -’jv  iiauw*  •_• 


Fart  JY > Conclusions 

236.  The  as*  theory  of  the  boundary  layer  for  turbulent  flow 
(inferences  (»)  and  (q))  is  governed  by* 

a,  The  continuity  of  decay  for  both  heat  and  momfentum 
effects  of  turbulence  froa  the  free-stream  * ) the  wall. 


b,  A better  compliance  with  the  requirement a of 
continuity  for  mass,  momentum  end  energy  in  regard  to  both 
the  mean  flow  and  the  turbulent  motion.  This  led  to  the 
discovery  of  a very  satisfactory  fora  for  the  function* 
p » Hr/h*  0Tsr  entire  boundary  layer. 

e,  Interpreting  the  experiments  of  f lev  over  heated 
bodies  ir*  tha  sens*  of  assigning  to  the  unknown  function 
m.  * A}  f Ar  * ovsr  the  boundary  layer,  the  border  values 

* l -w.  V *■  , and  estimating  ““aver.  for  various  flow 

W UAJLA  A bAwUC  • 


1 - The  use  of  Pr  ' = mPr  as  one  of  the 

essential  parameters  in  dealing  with  boat  transfer  problems  of 
turbulent  flow. 

e.  A distinction  between  two  characteristic  layer-widths 
(differing  for  pr  ! *9=  l 1)  which  are*  % for  viscous  effects, 
'Vt  for  heat  flux  density  . Sach  one  marks  a limit 

distance  from  the  wall  within  which  the  relevant  turbulent 
effect  ( ti  or  ) does  not  exceed  its  molecular  counterpart 

< %.-  )• 

237.  Based  on  these  five  points,  the  calculations  of  the 
temperature  distribution  and  heat  transfer  have  been  extended 
over  s range  ofs 


High  ysl  >uoa  O?  vuv  ?r -number  (References 


up  to  Pr  * =•  10-3 , 


b.  Low  Pi -numbers  (References  (q  ) , (m  ) , and  (a))  down 
tu  Pr'  = 10-3, 0) . * ' 

236.  The  former  category  corresponds  to  nonmetallic  liquids 
like  water  and  oils.  The  latter  covers  conditions  of  metals 
In  liquid  state  (mercury,  solten  laetals  and  alloys)  • The 
lower  range  becaioe  of  particular  importance  because  of  the 
necessity  of  removing  heat  in  large  amounts  from  nuclear 
furnaces,  at  a high  temperature  lave?,  by  applying  molten 
metals  as  coolants. 
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239.  The  shape  of  the  temperature  profile  appears  to  be 
highly  dependent  upon  the  generalised  Pr-numbe.  (?r*  * mPr).. 
whereas  it  seems  leas  affected  by  the  Re -number 1 However,  th« 
dimensionless  heat  flux  density  ( <•  / ‘t  ) is  found  more 
dependent  on  the  Re -number. 

240.  Almost  irrespective  of  both  parameters  (Pr»  and  Re)  the 
average  mean  temperature  for  the  turbulent  flow  in  a circular 
tube  sppearea  a z a distance  from  the  wall  equal  to  0.3  of  tha 
tube  radius  (References  (p)  and  (q)). 

2/4J..  As  for  the  task  laying  ahead  of  us  in  this  field,  let 
the  following  be  hinted.  A considerable  need  seems  to  subsist 
fer  exploring  still  more  by  systematic  experiments  the  profiles 
of  both  velocity  and  temper® ture  for  condition#  of 

turbulent  flow.  Increased  accuracy  of  results  might  allow  one 
to  reach  by^ inf erenoc  the  form  of  the  characteristic  function* 
m = ''V  / Wf  . This  would  be  an  important  contribution 
toward  a better  understanding  of  what  appears  to  be  a composite 
structure  of  turbulence. 


« 


62 


c>*l 


NAVOHD  Report  2956 


REFERENCES  * 


(a)  Reynolds,  0-, 


» A„  X.».  . 

vu  wiiO 


Extent  and  Action  of  tho  Heating  Surface 
Staam  Boiler8*.  Fr°,c*  Manchester  Lit.  and 

Phi  J.  . SOC.  VOI.  l4.  l£“ll  r\r»  "7  ..1  O J1  4.. 

Conectcd  Papar^.  Cambrldg^  EngTaAd  KIT  lV 

No.  lq.,  1900,  pp.  81-85. 


<b)  Stanton,  T.  E., 


(*}  rt0n  the  Passage  of  Heat  between  Metal 
Surficis  jnd  Li(juiu$  in  Contact  aith 

thss!?"  oShi1,  Trans*  R°y*  Soc.  A 190,  1897, 

PP . 07 '03. 


(ft)  "The  Relation  betwoen  Friction  and  Heat 

Transmission  in  the  Case  of  the  Motion  of  a 
Fluiu  over  a Surface  having  Temoerature 
differing  from  that  of  theFluid. " Diet. 
Appl.  Physics.  Vol.  1^  Art.  Friction.  VI  f 
40-4-1,  1922,  pp.  399-404. 

(c)  Prandil,  L., 

Relation  between  Heat  Exchange  and  Drag" 
Sine  Be z is hung  swischen  Wirmeaustausch  und 
Stromungswider  stand  der  Fliissigkeiten." 
Phys.  Zeitschr.  11,  1910,  pp.  IO72-IO78, 

((*)  Remarks  on  tha  Heat  Transfar  ir.  the  Tube." 
oemerkungen  uber  den  Wiirmeubergang  in 
Rohr."  Phys.  2,  29,  1928,  pp.  487-489. 

(d)  Taylor,  O.  X,, 


(*) 

<(>> 


"Eddy  Motion  in  tha  Atmosphere."  Phil. 
Trans.  Roy.  Soe.  A 215,  1915,  pp.  1-33. 


"The  Transport  of  Vorticity  and  Heat 
Fluids  in  Turbulent  Motion."  Proc. 
Soc.  A 135,  1932,  pp.  685—705. 


throu  gh 
Roy. 


* All  reference  titles 

• -VI  4 A-  * — A-  ..  . _ 

wawm  XU  «SL  L\J  UliO i* 

r&nslation. 


are  given  in  English.  • However,  for  \ 
languwaM.  the  original  title  follows  Its 


NAVORD  Report  2956 


REFERENCE?.  (CCNT?D) 


(•)  Stodola,  A., 

"To  the  Theory  of  Hast  Transfer-  from  Gases  or 
Liquids  go  Solid  Wails"  — 8Zur  Theorie  des 
Warmeubergangea  von  Geaen  Oder  Flussigkeiten 
an  rest*  -and «.n  sehweis.  Baustg.  1926, 

pp. 

243-  244.  1927,  pp.  193-196. 
it)  vcn  X=rman,  Th., 

(*0  "Turbulence  and  Skin  Friction."  Joum. 
IAS  I934#  PP-  1-20. 

(ft)  "The  Analogy  between  Fluid  Friction  and 
Heat  Transfer"  ASME,  Nov.  1939,  pn,  705- 
710. 

(g)  Elias,  F., 

"Ee&t  Transfer Kat  a Heated  Plate  to  Flowing 
Air"  — "Die  Warraeilbertragung  einer  geheiaten 
Platte  an  atr-<Jmcnda  Luft."  Auhandl.  Aachen. 

J**  WO.  Also  ZAMM,  i.929,  pp.  434-453, 

II  I93O,  PP*  1-14. 


(h) f  Colburn,  A.  P., 

"A  Method  of  Correlating  Forced  Convection  Heat 
Transfer  Data,  and  a Comparison  with  Fluid 
Friction."  Trans.  Am.  Inst.  Cham.  2ng.  Vol 
XXIX,  1933,  PP.  174-210. 

(i)  Lorens,  H.,  and  Friedrichs. 

"Contribution  to  the  Problem  of  Heat  Transfer" 
"Beitrag  t'un  Problem  des  Wlsrnieuberganges." 
Zeitachr.  teohn.  PVnr*ilr.  lQili  i dA_i Ao 

201-206,  376-377.  * ~ ' ' 

( J ) Kaye , W , A. , and  Furnas . C . C . . 

"Heat  Transfer  involving  Turbulent  Fluids", 
Induafcr,  Engin.  (Jhem.  Vol.  26.  193a.  dd.  78^- 
r 66. 


64 


it  A VO  HD  Report  2958 


RSraRKflCES  ( OCN? ! D ' 


(k)  Sioder . E.  N.,  and  Tate,  3.  E. , 


JSS!^Tr^!£*r  anlPra”ul‘,!>  D**°P  of  Liquids  in 

W^'il)i9"^35.“  ‘ “**  ‘ Vo1*  a8»  W36. 


(1)  isokert.  B., 


(m)  Reichardt,  H., 


( > oea&  Transfer  through  Turbulent  Friction 
Layer.."  NACA,  TSf  lQfcy,  Deo.  19^3,  53 

*£? ‘ 297!3288A*ti0n  fr0m  ZASM*  D*°*  19^°* 

(pO  "The  Influenoe  of  the  Flow  Near  the  W*n 
upon  sns  Turbulent  Heat  Transfer"  — "Dsr 
~r  ~nr  3 8 -®r  yandnahen  St-romung  auf  den 

tUrbulAnt^^  V2vm«  tv  ||  • » 1.  - __ 

a Si*  Mltteil.  Nr. 

1 of  Max  Planok  Inst.,  Odttingen  1950. 

63  pages.  * 


(i) 


(n)  von  Mi aea,  r.. 


Complete  Representation  of  the  Turbulent 

°£  Y«locity  in  Smooth  Duote"  - 
» t e ! lun g der  turbulenten 

£•«»  _. 

» * -/,>.(  two-cjLO. 

The  rounaatione  of  the  Turbulent  Heat 

r 01  2 °fundlagen  des  turbulente: 
WaraeUberganpa."  Arch.  gea.  Wiirmetechnik, 

VUA  v * inuiiar  1 C\  •?  a J.  9 

* 0*“w  - 7 a yy  e ; 


JL®2™  Ramarka  on  the  Laws  of  Turbulent  Motion  in 
pp!  317-327?  <tnain  Annlv*ra“T  Folume,  I9L1. 


(o)  Hauaen,  2. , 


"Representation  of 
Crenaralized  Power 


the  Heat  Transfer  in  Tubas  by 

Relations"  - "Daratellunr  de* 

+**  - — 


a wf 

vs 


NAVCKD  Report  2953 


ppygu  u-r~~  tcQ  ( P Q r^r  • TJ  ) 


: .Tmirn\, 

* f i » i 


W^raeiibergangee  in  Rohran  durch  versilg*5ieinarte 
Potenzbeziefcungan. " VDI  Verfahrenaleohn.iic,  1943, 
Nr.  4,  pp.  91-98. 


Li  sou* s ion  by, 

{*<)  Tress.  J La.  Inst.  Chem.  lag.  1940#  PP*  657- 


oo^. 


t A \ m M «.  • 

153*/  14  04*0.  ; 


69.  1947,  p.  958, 


{ q ) Siser , K. , 


(*4)  ’‘Temperature  fields  due  to  Viscosity  in 

Turbulent  Boundary  Layers*  *>--  "Reibungstea- 
peraturfelder  in  turbulenten  Grenzaohiohten. * 
Hitt all.  Hr.  8,  last.  Thermodyn.  and  Design 
Comb.  Motors,  STfi.  Zurich,  1949,  88  pages. 

(A)  "The  Turbulent  Heat  Transfer  in  the  Tube  at 
r very  Small  Frandtl  Humbers"  — "Der  turbu- 
lent* Warmeubergaag  im  Rohr  bei  sehr  Jclelnen 
Prandtl-Zahlen. " Aroh.  gas.  Warmeteohn. 

Nr.  10,  1951,  pp.  206-211. 


(r)  Koch.  R. , 


"Turbulent  Heat  Interchange  In  the  Tube"  »** 
"Turbulenter  Waroeausteusoh  io  Rohr."  Arch, 
ges.  Warmeteohn.,  Jan.  1950,  pp.  2-8. 


(s)  Lyon,  R.  N. , 


"Liquid  Metal  neat  Transfer  Coefficients."  Chem. 
Sag.  Progress,  February  1951,  PP*  75-79= 


(t)  Theodoridea,  Ph.  J. , 


"Turbulence  and  Heat  Transfer  at  Smooth  Walla 
for  a Flow  in  a Prismatic  Duot  or  along  a Flat 
Plate"  {Unpublished  Survey  of  Soma  Recent 
Developments),  1952,  72  pages 


66 


I 


NAVORD  Report  29.56 


ILLUSTRATIONS 


Figure  i.  Analogous  coefficients  { ) for  time  rate  of 

molecular  and  turbulent  transport  of  momentum, 
heat,  and  mass,  respectively  ...  .... 


2.  Ratio  of  turbulent  to  viscous  friction, 
relative  velocity  and  Quantity 

(at  different  Frandti  numbers)  in  terms  of 
nondimensional  distance  from  wall.  (Graph 
reproduced  from:  H.  Reichardt,  Mitteilung 

No.  3 of  Max  Planck  Institute.  GAtti™ger># 
1950,  Abb.  1)  . . 


The  integrand  function:  ^ ( ),  aiid  the  ratio  of 

turbulent  part  to  total  shear  stress  over 
nondimenaional  distance  from  wall.  (Graph 
reproduced  from:  H.  Reichardt.  ZAMM „ July  1951, 

Bild  2)  ...  ... 


Relative  velocity  ( * u/u  ),  and  cross-derivatives 
thereof  ( . over  a rectangular  channel. 
(Qraph  reproduced  from:  H.  Reichardt,  ZAMfc-  July 

1951.  Biid  1*.) ........  . . : . 


5.  The  function  A/f>*lx  and  the  additional  relative 
velocity  ( '•*-/  I At*  ),  over  the  whole  cross  section 

for  flow  with  pressure  drop.  (Graph  reproduced 
from:  li.  Reichardt,  ZAMM,  July  19*51,  Bild  5) 


6.  Comparative  plots  of  relative  velocity  with  and 
without  corrective  tom,  ( Graph  reproduced  from; 
E.  Reichardt,  ZAHM#  July  1951.  Bild  7) 


Complete  velocity  profile  of  _ o » in  a cylindrical 
tube  for  different  Re-numherg , (Graphs  reproduced 
from:  H.  Reichardt,  Hittsilung  :o.  3 of  Max  Planck 

Institute.  Gottingen,  195- , Abb.  3 and  Abb»  ij.) 


8,  Distribution  of  dimensionless  heat  flux  density, 
aa  well  aa  shear  stress  , in  tube  ever  the 

distant:*  from  wall  (for  Re  = (3)  104).  (Graph 

reproduced  from  II , Reichardt;  Grundlagsn  do s 


turbulonten  Warrceuberganeres . Arch.  Gas.  Warmetechnik 

1951,  Heft  6/7,  Fig.  5)  . - 
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Figure  9* 


10. 


11. 
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Temperature  profiles  ower  nondimansional  velocity 
('Mt  and  ov«r  nondimeneional  distance  fro  in  wall. 
(Graphs  reproduced  from  H*.  Reiohardti  iii  tteilung 
Ho.  1 of  Has  Planck  Institute.  Gottingen,  1950, 

«*-u  ‘ C.  a ivv  •*  \ 

AUOi  V|  *****  »mm*  ; / 

Temperature  distribution  (£“)  vsrauB  7eioeifc7  yf). 
Turbulent  flow  In  tubes  (aecond  approximation)  for 
Re  = C s ) 10H*  (^=  = SOO)  at  Tsuious  ?randtl-mu»b»p*. 
(Cfcatph  reproduced  from  H.  Reichardt*  Orund lagan 
d««  turbulsnten  Wairaeuher  ranges,  Aroh.  oea. 
~£raei.chnik  19>1*  Heft  6/7#  Fig.  3)  - 


ft- 

Tcssporatur©  distribution  (-J' ) to r bus  dimensionless 
distance  from  wall  { i/X  ) . 


n — e i.  p — . x / .«  ...J 
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He  ■ (3)  10^-3  at  various  Pr and tl -numbers . ( Graph 
reproduced  from  E.  Reichardt t Grundlagcn  das 
tupRulenten  tflrmettbsrgangee,  Aroh.  Qee.  Wirmsteoha. 

1951,  H.  6/7,  Pig.  i*)  . . 
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rigur®  x is  to  illustrate  the  accepted  ways  or  expressing  the  effects  of  turbu- 
lent agitation  on  the  transport  of  nmaontun-  enthalpy  and  mass.  It  Is  in  close 
analogy  to  the  formalism  inspired  by  Newton  aiii  Fourier  for  viscous  and  heat 
conducting  effects  of  thermal  agitation.  In  addition.  Figure  1 includes  the 
usual  breaVdoryn  of  fluctuating  physical  quantities. 
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Ratio  of  turbulent  to  viscous  friction,  rela 

tive  velocity  ( ) , and  quantity  a,^  (at  differ 
ent  Prandtl  numbers)  in  terms  of  the  non-dimen- 
sional distance  from  the  vail. 

(1) 

(Graph  reproduced  from:  K.  Reichardt,  Mitteilung 

No.  3 of  Max  Planck  Institute,  C-ottingen,  195° , 

Abo.  l) 
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shear  stress  ever  non-dimensional  uis 
tance  ( h ) from  vail. 


/ ^ V 
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(Graph  reproduced  from;  H . Beichardt 
ZM£4,  July  1951,  Bild  2) 
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Af 

The  function:  //*’?*  , and  the  addi- 
tional relative  velocity  ( ■ui/n*'  ) 

over  the  whole  cross-section  for  flow 
with  pressure  drop. 


(4) 

(Greph  reproduced  from:  H.  Beichardt, 

July  1951.  Bild  5) 

— ~ — * ^9 


Figure  5 


Comparative  plots  of  relative  velocity 
vith  and  vithout  corrective  term „ 

(5) 

(Graph  reproduce,  from:  H.  Beichardt , 

ZAMM,  July  1951,  Slid  7) 


Figure  6 
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Distribution  of  dimensionless  heat  flux 
density  %/  c^. c as  well  as  shear  stress 
*Ao  1**  tube  over  the  distance  from 
wall  (for  Re  * (3)  104  ) 

(Graph  reproduced  from  H*  Reichardt : 
Grundlagen  des  turbulenten  WUrmettber- 
ganges,  Arch.  ges.  Wfirmetechn.  1951,  H. 
6/7,  Pig.  5,  p.  136) 


Figura  8 
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Temperature  profiles  over  non-dimensional  velocity 
( ) and  over  non-dimensional  distance  from  vail. 

(7) 

(Graphs  reproduced  from?  H.  P.ei=b±.rdt,  Mitteiiung 

JMO-  Of*  H* T Plo«/»lr  * — IrtCA 

* ^ 1 w*  vMVV  j WV%W  j J*J  j 

Abb.  6 and  Abb.  7) 
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Temperature  distribution  (v*  ) versus 
velocity  ( r ) • 

Turbulent  flow  in  tube  (second  ap- 
roximation)  for  Re=(3)l0^  = BOO) 

t various  Prandtl-numbers • K 

(Graph  reproduced  from  H.  Reicnardt 
Grundlagen  des  turbulenten  WarmeUber- 
ganges,  Arch,  ges.  Wftrmetechn.  1951, 

H.  6/7,  Fig*  3-  p.  135) 
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0.2  0.4  0.6  y 0.8 


Temperature  distribution  (v)  versus 
dimensionless  distance  from  wall  (y/r). 

Turbulent  flow  in  tube  (second  ap- 
proximation) for  Re=( 3)10  at  various 
Prandtl-numbers * 

(Graph  reproduced  from  H.  Reichardt: 
Grundlagen  des  turbulenten  Wftrmeiiber- 

granges « Arch*  gesj  Wansetechn*  1951*  H« 
6/?:  FIs  4*  p.  136) 
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